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ABSTRACT 

We study the statistical mechanics of random surfaces generated by iV x iV one- 
matrix integrals over anti-commuting variables. These Grassmann-valued matrix 
models are shown to be equivalent to N x N unitary versions of generalized Penner 
matrix models. We explicitly solve for the combinatorics of 't Hooft diagrams of 
the matrix integral and develop an orthogonal polynomial formulation of the sta- 
tistical theory. An examination of the large N and double scaling limits of the 
theory shows that the genus expansion is a Borel summable alternating series which 
otherwise coincides with two-dimensional quantum gravity in the continuum limit. 
We demonstrate that the partition functions of these matrix models belong to the 
relativistic Toda chain integrable hierarchy. The corresponding string equations 
and Virasoro constraints are derived and used to analyse the generalized KdV flow 
structure of the continuum limit. 
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1 Introduction and Summary 



Matrix models |T| provide a useful framework in which to extract quantitative informa- 
tion about complex statistical systems given only their general symmetry properties. The 
key feature which enables this description is the existence of universal behaviour in these 
models. Over the years random matrix models have been utilized in a wide range of stud- 
ies in condensed matter systems ||, quantum chromo dynamics ||, and low- dimensional 
string theory [Q, || . They are also useful for solving various statistical mechanical models 
and combinatorial problems ||. One interpretation of random matrix ensembles is in 
terms of generators of 't Hooft diagram expansions which correspond to discretizations of 
compact Riemann surfaces. The continuum limits of such discretizations are realized as 
phase transitions in the statistical mechanics of the matrix model and characterized by 
universal data at the critical points. 

Despite the success of matrix model technology over the years, there remain several 
issues in random surface theory which are not adequately addressed by standard ap- 
proaches. The simplest, and one of the most fundamental issues regards a self-consistent 
definition of the vacuum of pure two-dimensional quantum gravity. The conventional 
matrix model methods for investigating this theory lead to inconsistencies. For instance, 
bosonic matrix integrals are typically ill-defined in the physically interesting region of 
parameter space corresponding to the non-perturbative continuum limit. Furthermore, 
they are unable to describe the coupling of gravity to conformal matter fields of central 
charge c > 1. It has been suggested J?| that above this c = 1 barrier the two-dimensional 
geometry degenerates into a tree-like or branched polymer phase. 

In this paper we will study a model of discretized random surfaces which is in part 
motivated by such issues. It is parameterized by a zero-dimensional action involving 
matrices with Grassmann- valued elements |§. This model is known as the adjoint fermion 
one-matrix model and it is defined by the partition function 

Z N = J dip dip e NtlV ^ (1.1) 

Gr(AT) c 

where 

n#) = E? {^) k (i.2) 
fc=i k 

is a polynomial potential of order K < oo, and dtp dip = flij dipij dipij is the standard 
Berezin measure on the complexified Grassmann algebra Gr(iV) c = Gr(iV) ®r C. For 
any matrix X the trace is normalized as trX = J2iXu. The matrix elements ipij and 
ipij, i,j = l,...,N, are independent, complex-valued Grassmann numbers with the usual 
rules for complex conjugation, ip*j = —ipij, ip*j = —ipij, and {ipijipki)* = iph^ij (This 
ensures that the generating function (|1.1|) is real- valued) . The large N limit of this 
model has been studied in 0-0 and shown to exhibit the range of critical behaviour 
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seen in the usual bosonic one-matrix models. The topological expansion of the matrix 
integral ( [Ll|) was considered in [10] and some novel critical behaviour was observed in 
[|TT| . The fermionic matrix model can also be coupled to an ordinary complex bosonic one 
to produce a supersymmetric matrix model, which has been used to describe solutions of 
certain combinatorial problems [12| and also to generate statistical models of branched 
polymers |13 . 



The most attractive feature of the fermionic matrix model (|1.1|) is that the integration 
over Grassmann variables is always well-defined and the theory yields finite, computable 
observables. The dimension N of the Grassmann matrices acts as a cutoff on the number of 
terms in the partition function (|1 . 1| ) which counts all fat-graphs generated. The fact that 
the partition function is a polynomial in the coupling constants for finite N, coupled 
with the fact that inner products in the Grassmann integration measure are effectively 
calculable, leads to explicit expressions for statistical quantities in this class of models. In 
the following we will use this property to solve analytically the combinatorial problem of 
counting the dynamical triangulations generated by the fermionic partition function via 



an explicit determinant representation of (1.1 



Since the generating function fll.l]) is always a well-defined convergent object at finite 
N, one might expect that its large N limit is also well-defined. In perturbation theory, the 
Feynman rules for fermionic matrices include a factor of — 1 for each fermion loop thereby 
resulting in cancellations between large numbers of Feynman diagrams. This observation 
has led to the conjecture that the genus expansion of the matrix integral (|1 . 1| ) in the 
large N limit is an alternating series which may be Borel summable. The evidence for 
this striking property of the adjoint fermion one-matrix model is its intimate relationship 
with generalized Penner matrix models @, [10], || which are defined by the Hermitian 



matrix integrals 

Z* a = J d(j)e Ntr(vW +a lo g< p) (L3) 
u(N) 

where d(f) = Yiij dipij is the Haar measure on the Lie algebra u(N) of the N x N unitary 



group U(N). It has been argued that the matrix model (|1.3|) is equivalent to ( |1 . 1|) for 
a = —2, with the Hermitian matrix <fi identified with ipip. This has been established in 
the large iV limit by a saddle-point computation on ~ ipip 0, and by showing that their 
Schwinger-Dyson (loop) equations are identical order by order in the large N expansions 
of the matrix integrals || |Hj . However, for a < the integration over Hermitian matrices 
in ( p.. 3D is not well-defined because of the logarithmic divergence at <fi = 0. In this region 
the integral can only be defined by analytical continuation and the matrix model only 
makes sense at N = oo. This continuation produces complex- valued endpoints for the 
support of the corresponding distribution of </> eigenvalues, analogous to what naturally 
occurs in fermionic matrix models || It is this complexification of the support of the 
distribution of eigenvalues that is asserted to result in an alternating genus expansion. 
Moreover, it can be argued [10] that the topological expansion coincides, modulo these sign 
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factors, with the usual Painleve expansions of Hermitian matrix integrals with polynomial 
potentials [|J. 

The near equivalence of the genus expansions strongly suggests a correspondence be- 
tween the Feynman graph expansions in the bosonic and fermionic models. One way of 
understanding these issues is to appeal to a simplified version of the matrix integral ( |1 . 1| ) , 
namely a fermionic vector model defined by anticommuting vector elements ipi and 
i = 1, . . . , N JT4|. A fermionic vector model is related to a bosonic vector model, defined 



with the same polynomial potential V, by a simple mapping of its Feynman diagrams and 
an analytical continuation JV m — | of the vector dimension which results in a bosonic 
perturbation series with a cutoff on the number of terms generated. This analytical con- 
tinuation can be understood in terms of a complex-valued saddle-point of the bosonic 
model in the large N limit, leading to an alternating "genus expansion" which produces 
a well-defined, unique function WM. As bosonic vector models generate random models 



of branched polymers, the fermionic vector model thereby produces a statistical theory of 
branched polymers whose continuum characteristics are equivalent order by order in large 
N to the conventional models, but whose double-scaling expansion is a non-perturbatively 
well-defined function. 

The arguments given in [1(J are based on the class of models with odd polynomial 
potential, i.e. V(—<f)) = —V(<f)). In these cases, the endpoints of the support contour 
of the spectral density are located symmetrically about the origin on the imaginary axis 
in the complex plane || [|. Generic polynomial potentials are difficult to treat because 
the corresponding support contours are located asymmetrically in the complex plane and 
the loop equations become cumbersome to deal with. Moreover, although the equations 



of motion for the matrix models and (|l.q ) are identical at each order of the large 
N expansion, beyond the leading order they must be solved with different boundary 
conditions and the non-perturbative solutions are not the same in the two cases JED] . One 
of the main results of this paper will be the clarification of the role played by generalized 
Penner models in fermionic matrix integrals. The analytical expressions for the partition 
function ( |1.1| ) at finite N will be interpreted in terms of Toeplitz determinants, which 
will naturally lead to a proof of the equivalence of the fermionic one-matrix model and a 
unitary version of the Penner one-matrix model, 

Z%= J [dU}e Ntr ( v ^- l ^ u ) (1.4) 

U(N) 

where [dU] is the invariant Haar measure on the unitary group U(N). The important 
feature of this equivalence is that it holds for any matrix dimension N. The finiteness 
of the fermionic matrix model at finite iV is captured by the compactness of the field 
variable in ( |1.4|) . As we will see, the logarithmic term in the effective potential of ( |1.4|) 
gives certain restrictions on the configurations which contribute to the partition function 
and its observables when the unitary matrix model is treated from a group theoretic 



4 



perspective. These restrictions naturally reproduce the basic characteristics of the matrix 
integral (|1.1| ), and moreover show precisely how the analytic, functional forms of the 
models (|1 . 1| ) and ( |1.4| ) are equivalent for any value of N. We will also prove directly that 



the loop equations of the fermionic and unitary matrix models are identical, and that 
they admit the same, perturbative Gaussian boundary conditions. It is a remarkable fact 
that a unitary matrix model such as (|1.4j) admits an interpretation in terms of random 
surfaces. 

The biggest advantage of the mapping of the fermionic matrix model to a unitary one 
is that we now have an eigenvalue model to analyse. Grassmann-valued matrices are not 
diagonalizable, and a direct analysis of fermionic matrix models is only possible using 
the method of loop equations || [|. In this paper we shall exploit this eigenvalue rep- 
resentation to give a systematic and detailed account of the exotic properties possessed 
by fermionic matrix models. In particular we will develop a generalization of the orthog- 
onal polynomial technique on the circle which is necessary to study the unitary matrix 
model defined by ( |1.4j) , and describe some of their functional analytic properties. These 
methods will naturally lead, irrespective of the parity of the potential V, to a precise 
investigation of the double-scaling limit of the fermionic matrix model which determines 
a non-perturbative, all-genus continuum theory of two-dimensional quantum gravity cou- 
pled to matter. We will show that the fermionic matrix model leads to a Borel summable 
genus expansion, in contradistinction to Hermitian one-matrix models. Unlike the vector 
models, the diagrammatic, finite N situation is not quite so simple in the case of fermionic 
matrices, but the final qualitative result will be the same. In particular, we shall obtain 
a non-perturbative definition of pure two-dimensional quantum gravity which is equiva- 
lent to the usual models order by order in the genus expansion, but which is given by a 
well-defined unique generating function. This non-perturbative model is what we shall 
call "fermionic quantum gravity" in the following. 

Another advantage of the eigenvalue representation is that it will allow us describe 



the fermionic one-matrix model ( |1 . 1| ) as an integrable system. Generally, matrix models, 
being examples of exactly solvable systems, are intimately related to certain reductions 
of well-known integrable hierarchies of differential equations jl5fl . In the usual bosonic 



Hermitian one-matrix models, describing two-dimensional quantum gravity, the partition 
functions are related to the r-functions of the integrable Toda chain hierarchy. In the 
following we will see that the fermionic gravity theory is related to the r-function of a 
particular deformation of this hierarchy known as the relativistic Toda chain. We can 
therefore understand many of the novel aspects of the fermionic theory, such as what 
sort of matter coupling it involves, in terms of this deformation. Furthermore, we will 
develop an operator theoretic approach to the equations of motion of the matrix integral 
(|1 . 1|) and gain in this setting a precise picture of how this deformation leads to a well- 
defined Borel summable partition function of the statistical model. This sets the stage for 
a description of the integrable differential hierarchies satisfied by the continuum theory 
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which generalizes the usual KdV flow structure of gravity j|, H . It will also indicate that 
the partition function (|1.1|) in the continuum limit serves as a concrete, non-perturbative 
definition of two-dimensional quantum gravity. 

The organization of the remainder of this paper is as follows. In section 2 we will 
derive an analytical expression for the partition function ( |1.1| ) and use it to describe the 
solution to the problem of counting fermionic fat-graphs. We shall see that the matrix 
integral generates a novel class of discretizations of Riemann surfaces whose polygons are 
two-colourable. We will also briefly describe what sort of matter-coupled gravity theory is 
represented by the fermionic one-matrix model, and prove its equivalence to the unitary 
one-matrix model (|1.4j) . In section 3 we will present the formal orthogonal polynomial 
solution of the fermionic matrix model. In section 4 we will describe the large N limit of 
the model from the perspective of these orthogonal polynomials. We will find that the 
adjoint fermion one-matrix model actually possesses a sort of "internal" branched polymer 
critical behaviour which can be understood in terms of the fermionic characteristics of 
its fat-graph combinatorics. We shall also demonstrate that the usual universality classes 
of two-dimensional gravity arise, and that they cannot be smoothly connected to the 
branched polymer phases of the theory. Of course, the polymer effect is washed out 
at iV = oo because the surface effect is of higher order in N. We also construct the 
genus expansion of the theory and establish its Borel summability explicitly. In section 5 
we present the operator approach to the fermionic matrix model and use it to interpret 
the partition function as a r-function. We derive the Virasoro constraints satisfied by 
this r-function, which at the same time establishes explicitly the equivalence of the loop 
equations of the matrix models (|1.1| ), ( p..3|) with a = —2, and ( |1.4j ), in the large N 
expansion. We then use these descriptions to give an functional analytic explanation for 
the nature of the topological expansion in the fermionic one-matrix model. Finally, in 
appendix A we describe properties of observables of the statistical model ( |1.1| ) using the 
unitary matrix model and its orthogonal polynomials, while in appendix B we show that 
the orthogonal polynomials of the Gaussian fermionic one-matrix model are given in terms 
of confluent hypergeometric polynomials. 



2 Determinant Representations of the Partition Func- 
tion 

In this section we will present an analytic solution to the problem of counting fermionic fat- 
graphs. We will obtain a closed expression for the generating function of such discretized 
random surfaces, and at the same time we will be naturally led to the equivalence of the 
fermionic matrix model and the unitary matrix model which will occupy the bulk of our 
analysis in subsequent sections. We will also obtain in this way a geometric picture of 
what sort of theory of random surfaces is represented by the adjoint fermion one- matrix 
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model and also of the precise role of the corresponding Penner interaction. 



2.1 Hermitian Representation 

One intriguing feature of the fermionic one-matrix model ( |1 . 1|) is the extent to which it 
can be solved at finite- N. This is a consequence of the convergence of the integration over 
Grassmann variables, in contrast to the case of bosonic matrix models. The combinatorics 
of the fat-graph expansion of the fermionic partition function can be determined by 
using the observation of [JTTJ that the partition function and observables of the matrix 



model ( |1 . 1|) can be evaluated explicitly by introducing two N x N Hermitian matrices 
and A defined by 

dX 



^Mn-rr) . Mo)-- j _^_ e ^ (21) 



u(N) u(N) 



By inserting these definitions into (|1 . 1|) , the fermionic partition function can be written 

as 

dX 



e iVtrV(0) j (lc(ll , j c *tr 

1 



(27T) 



N 2 



u(N) Gr(A r ) c u(N) 

' d<pdXe NtvV ^ +itlX ^ det N {-iX) 

u{N) 

d<f>e N ^det N (-^)5(4>) (2.2) 

u(N) 



so that 



to* (£) e"-(« |^ (2.3) 
A similar expression can also be derived for the correlators of the fermionic one-matrix 



model |TT|, |9| . However, it is difficult to write down an explicit expression for ( |2.3| ) which 
is informative and amenable to analysis. 

To evaluate the determinant in ( j2.3| ) explicitly, we shall instead write the fermionic 
matrix model as a random theory of the two Hermitian matrices <fi and A introduced above, 
keeping only the second line of the identity fl2.2|) . In this way, the fermionic one- matrix 



model can be written as the Hermitian two-matrix model 

N N(N+1) 



(27T) 



N 2 



d\ e #tr(v(0)+log(-iA)+iA<£) (2.4) 



u(N) 



where we have rescaled A^ — > N ■ A^-. The doubling of degrees of freedom is required to 
compensate for both the matrix ip and its adjoint ip. The Hermitian two-matrix model 



(2.4) describes two-dimensional discretized gravity with a single type of matter state at 



each of the vertices. It does not admit the usual Z 2 Ising symmetry of a classical spin 
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lattice that characterizes the unitary minimal models of rational conformal field theory 
[[| . Thus the matter states in the worldsheet interpretation of the fermionic matrix model 
are more complicated degrees of freedom than just simple Ising spins or other conventional 
types of conformal matter. The logarithmic nature of the effective two-matrix potential 



in (|2.4j ) is just another reflection of the connection with Penner matrix models. The 
original Penner model |16[ (gk 



for k > 1, g\ = —1 and a = 1 in ( |1.3| )) we used 
to calculate the virtual Euler characteristics of the moduli spaces of compact Riemann 
surfaces. For a > the generalized models are intimately related to the coupling of 
gravity to matter fields of central charge c = 1 JTOf- Note that the Penner interaction 
in ( |2.4| ) is well-defined. We may refer to the fermionic matrix model as a model of two- 
dimensional quantum gravity interacting with "Penner matter" . As we will see later on, it 
is this worldsheet model which leads to a well-defined nonperturbative theory of quantum 
gravity in two-dimensions which we will call "fermionic gravity" . 



The partition function (|2.4j ) can be written in terms of a double eigenvalue distri- 
bution by diagonalizing the Hermitian matrices <ft = U diag(0i, . . . , 4>n) U> and A = 
V diag(X\, . . . , Ajv) by unitary transformations, where <pi, A; G M are the eigenvalues 
of 4>, A. Computing the Jacobian for the change of integration measure and using the 



Itzykson-Zuber formula [[L7| to integrate out the unitary degrees of freedom, we arrive at 
the eigenvalue model 



N 00 

ZN = c N f[JJ dfa d\ X? e NV ^ +lNX ^' A(0 1; . . . , <M A(Ai, ...,X N ) 



(2.5) 



where c N = {-tN 2 ) N ^ N+l ^ 2 /(2n) N Un=i n\ and 



A(xi, . . . , x N ) = det 



J- 1 



i<j 



(2.6) 



is the Vandermonde determinant. Using (|2.6|) we can then write the partition function as 



Z N = c N N\ det [(f~ l , A^ 1 ) 
where we have introduced the inner product 

oo 

(F(0) , G(X)) =JJd<j)dXX N e NV ^ +mx * G(X) 



(2.7) 



(2.8) 



on the vector space C[0] ®C[A] of complex coefficient polynomials in (0, A). A remarkable 
feature of the fermionic matrix model is the extent to which its inner product ( |2.8|) can 
be computed. For arbitrary polynomial potential, integrating over A first in ( |2.8|) leads 
to the result 



N+k 



F(<j>)e 



NV{<j>) 



if>=0 



[2.9) 
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which is valid for any analytic function F((f>) and any integer k > 0. Using ( [2.9|) , the 
representation (|2.7|) of the partition function as the determinant of the moment matrix 
of the inner product (|2.8|) thereby yields 



•l)^ 2 det 

hi 







x N+j-i 




) 1 


<t>=0 



(2-10) 



where, for any integer K > 1, [N]k denotes the largest integer less than or equal to N/K. 
Here and in the following a prime denotes differentiation. The expression fl2.10| ) is the 
desired explicit expansion of the determinant operator in (|2.3|). 



2.2 Combinatorics of Fermionic Ribbon Graphs 



In contrast to the determinant form ( |2.3| ), the representation ( 2.1(J| ) leads to an explicit 
expression for the perturbation series of the fermionic one-matrix model. For a polynomial 
potential (|1.2| ) of order K, the determinant in (|2.10|) may be evaluated explicitly by 
expanding the exponential function in its Taylor series and applying the multinomial 
theorem to each term in the expansion to get 



d 1 



[Lh [L]k 

(%) £ E ••• E 



k 2 =0 



K 



9l 



■-Kk K )\k 2 \ 



x S vn^ { 9l y, 



(2.11) 



for any integer L > 0. Normalizing by the Gaussian model for which V(4>) = g\<p and 
^Gauss _ (_ij[N] 2 (jv^^^ 2 , the partition function ( |2.10| ) may be written as 



^ n N = det 
Z% auss i,j 



'[N+j-i] 2 [N+j-i] K 

E ■•• E 

k 2 =0 



(N + j 



K 

n 



9e 



fcf 



(N + j-i-Et tk t )\ k 2 \ ■ ■ ■ k K \ \£W-i ( 9l y 



(2- 



2) 

We can extract the sums in (|2.12 ) out of the determinant by using the multilinearity of 
the determinant as a function of its iV rows to get 



vGauss 
N 



E 



fcW..,fcW>0 



x det 

i,3 



E 

...,k 

(N + j-i)\ 



k£\...,k^>o 



K N 

n n 

m=2 n- 



K 



n© ([N+j 



Hi 



Si) 



(2.13) 



{N + j-i-Etikly. 1=2 

where denotes the step function with the convention 0(0) = 1. 

The expression ( 2.13Q leads to a relatively straightforward expansion of the partition 



function in powers of the coupling constants g#, £ = 2, . . . , K. It represents the formal 
solution to the problem of counting the (connected and disconnected) Feynman-'t Hooft 
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diagrams of the adjoint fermion one-matrix model. It is a closed formula, in terms of 
a sum over a large set of integers, for the generating function of ribbon graphs of the 
fermionic matrix model. A fermionic ribbon graph has much more structure to it than a 
conventional fat-graph, and so we shall now formulate the rules for generating them from 
the matrix integral. The fermionic propagator is (ipijipki) Gauss = 5u5jk/Ngi which we 
represent in the usual way as a double line, along with an orientation defined by an arrow 
pointing towards the if> matrix. Vertices are likewise given an orientation by assigning 
an outgoing arrow for a ip line and an incoming arrow for a if) line. Fat-graphs are now 
drawn with the rule that only if) and if> lines can contract. Each such graph is dual to 
some discretization of a Riemann surface in the standard way. The surface may be two- 
coloured by shading the triangles in a polygon if and only if their outer edge is dual to a 
propagator with an incoming arrow. In this way black edges are always associated with 
if> matrices and white ones with if). Since the fermionic propagator only connects ^'s with 
if>'s, it follows that every such discretization can be two-coloured. Thus the fat-graphs 
generated by the fermionic matrix integral can be obtained by drawing all discretizations 
associated with fermionic /c-point vertices in terms of the corresponding Hermitian 2k- 
point vertices, and keeping only those graphs which are two-colourable. The even parity 
of the Hermitian potential V((f> 2 ) is required because only discretizations with even-sided 
polygons can be two-coloured. 

However, the Grassmann nature of the generating matrices yields important changes to 
the rules for counting such triangulations. For a given vertex in the Wick expansion of the 
matrix integral, we fix a chosen if) line. Lines are joined into propagators with a left-handed 
orientation. Every time a if) line is joined into a if) line (rather than the other way around), 
we "twist" the ribbon and thereby reverse its orientation. To each such twisted line, we 
associate a factor of —1. This standard sign factor for fermionic fields yields significant 
reductions in the overall number of ribbon graphs which are actually generated by the 
fermionic matrix model, because of the cancellations which occur, for example, between 
fermionic diagrams that are topologically the same but which are twisted relative to one 
another. Thus topologically equivalent diagrams do not necessarily add up, but may 
have the effect of cancelling each other. As a simple example, consider the quadratic 
potential, K = 2. There are two graphs which contain only a single four-point vertex, 
but they contribute with equal magnitude and opposite sign to the generating function. 
This vanishing contribution is the coefficient of the g2f2N(gi) 2 term in the perturbative 
expansion of the partition function. Symbolically, it is given by the Wick expansion of 
the Gaussian correlator 

< tr <^> 2 L„„ = Q - CO = 

y^j (2.i4) 

Note the twist in the second propagator of the second graph in ( |2.i4j ), which induces the 
relative minus sign. Notice also that the usual toroidal four-point graph does not appear 
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in Q2.14Q , because it cannot be two-coloured. This vanishing contribution is also easily 



found from the analytical formula ( [2. 13] ). The leading term is 1 and it comes from the 
configuration whereby all /c's are 0. The terms of order g 2 come from the configuration 
whereby only one k is no n- vanishing and equal to unity. In the determinant in Q2.13|) , the 



only O's which are not equal to unity are those which appear in row i with fcM = 1. The 
determinant therefore vanishes, reproducing the graphical result ( [2.14| ). The remaining 



fat-graph combinatorics are now readily determined from the graphs of the Hermitian one- 
matrix model with even potential V(4> 2 ) by keeping only those which are two-colourable 
and incorporating the appropriate twists. The determinant formula (|2.13|) is readily seen 



to reproduce the correct combinatorics, with appropriate minus signs coming from the 
determinant^. 

These arguments readily generalize to arbitrary polynomial potentials. Generally, the 
total number of 2£-valent vertices in a given ribbon graph T as determined by the formula 



(|2.13|) is the power of gi which is given by 

N 



n*(r) = £*r (2.15) 

n=l 

and, because of Fermi statistics, this number is bounded as 

TV 

0<n 2£ {T) < £[2iV-n]< (2.16) 



71=1 



It follows that there are only a finite number of fermionic fat-graphs generated. The 
twisting of ribbons and also the finiteness of the perturbation series may be attributed 
to the inclusion of a logarithmic interaction —2 log</> in the corresponding Hermitian one- 
matrix model with potential V{<p 2 ). This role of the Penner potential will be demonstrated 
explicitly in the next subsection. 

The determinant expansion Q2.13 ) is therefore an analytic expression of the combina- 
torial formula 

Zn v (-i)^ (N 9l )*<r> x ( g e Y 2l{r) ( n 

z^r* r |Aut(r)i H W/ 1 • j 

where the sum runs through all (connected and disconnected) two-colourable fat-graphs V 
with n 2 ^(r) vertices of order 2£ bounded as in ( |2.16| ), and x(T) is the Euler characteristic 
of T. The alternating sign factor in (|2.17| ), with t(T) the total number of twisted ribbons 
of T, is due to Fermi statistics, while Aut(r) is the automorphism group of the unmarked 
graph T. The sum over fat-graphs in (|2.17|) is finite, since the maximum number of vertices 
that a given diagram can have is 

K N 

% ax = ^^[2iV-4 (2.18) 

£=2 n=l 



1 See [flUl for related considerations. 
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Each such ribbon graph is dual to a tessellation T* of a Riemann surface of Euler char- 
acteristic x(r) by ri2i(r) 2£-valent tiles for i = 2, . . . , K. Because of the twist factors, 
the overall combinatorial numbers generated by the fermionic matrix integral will be 
drastically reduced. Moreover, depending on the precise details of the potential V, the 
perturbative or topological expansions of the matrix model may be alternating series. The 
above considerations show that there is clearly no simple mapping between the Hermitian 
and fermionic matrix models. These features will have remarkable implications later on 
for the topological expansion of the matrix integral (|1 . 1|) . 



2.3 Unitary Representation 

We will now derive an alternative "dual" representation of the perturbative expansion of 
the fermionic partition function. For this, we note that the derivatives in (|2.9| ) can be 
evaluated in terms of a contour integration about the origin of the complex plane as 



(F(#) • A') = - (^) (N + k)l / dz [Z z Z +1 (2.19) 

2=0 



F(z) e NV & 



The partition function ( |2.10| ) may then be written as 



N / , • N N+k 



fc=l \ 

where we have defined the N x N matrix 



Z N = c N N\H(-(-j (N + k-l)\\ detA (2.20) 



r e NV{z) 

A H = f dz pvTT^TT ( 2 - 21 ) 

2 = 



The crucial observation now is that the matrix elements ( [2.21 ) depend only on the differ- 
ences j — i of row and column positions, i.e. A^- = A,_j. Matrices with such a symmetry 
are known as Toeplitz matrices and their determinants can be evaluated in terms of av- 
erages over the unitary group ]n|. The Toeplitz determinant det [A?-i] appearing in 
( |2.20| ) is associated with the Laurent series of a function A through the definition 



A{uo)= £ A n u n (2.22) 

n=— oo 

The matrix elements A n are then interpreted as the Fourier coefficients of A on the unit 
circle |u;| = 1, 

A n = cf-^- uj- n A(u) (2.23) 



Using the fact that a determinant is a linear function of each of its rows, it follows that 



at ^ 

det A = J] f A(uj k ) uo k k A(wi, ...,u N ) (2.24) 
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Replacing the right-hand side of (|2.24j) by its average over all permutations of the inte- 
gration variables u>k yields 



det A 



2 



1 r du 

mn/j^^K"' »»)' < 2 - 25 > 



which shows that the Toeplitz determinant is given by the unitary matrix integral 

detA= (2 ^" = (Li)/2 J [dU] det A(U) (2.26) 

U(N) 



Upon substituting ( 2.21|) into ( |2.22| ), we can interchange the sum and contour inte- 



gration provided we take into account that the latter vanishes for N + n < 0. We then 
have 

A(u) = i dz = 2ttz e Mv»-i°g-j (2.27) 

J uj" [z — uj) 

2 = 



From ( 2.20|) , (|2.26|) and Q2.27| ) it follows that the partition function of the adjoint fermion 



one-matrix model (|1 . 1|) is completely equivalent, for any dimension N, to the N x N 
unitary matrix model 

Z N = k N J [dU] (2 . 28) 

U{N) 

where k N = n n ^ n+1 ^ 2 /(2n) N< ^ N - 1 ^ 2 I\%=i( N + n- 1)!. We see that the fermionic one- 
matrix model for any iV is completely equivalent to both a Hermitian two-matrix model, 
and a unitary one-matrix model, both of which involve Penner-type interactions. In 
particular, the unitary representation clearly demonstrates that the equivalence with the 
Hermitian Penner one-matrix model is not true at finite N, but rather only at iV = oo. 
Note that the matrix integral ( [2.28] ) is perfectly well-defined at finite N, but it involves 
an integrand which is not a real-valued potential. 

With the partition function given as an integral over the unitary group, we can clarify 
the geometric role of the Penner-type potential that characterizes fermionic matrix models, 
and also give an alternative method for a perturbative expansion in terms of the coupling 
constants in the potential V. For this, we integrate over the U(l) factor of the unitary 
group U(N) = U(l) x SU(N)/Z N and leave an SU(N) integral, 

Z N = k N J d9 [dU ] e »(*ne ie u )-iNe) = ^ j d6 e"^ J [dU ] e NtrV ^ u ^ (2.29) 

U(N) SU(N) 

The role of the U(l) integration over 9 is to restrict the terms in the SU(N) integration to 
those involving only iV 2 factors of Uq. Equivalently, we may expand the invariant function 
which is the integrand of ( 2.29|) in U (N) characters Xft as 



e NtlV W=J2 c *Xn(U) , U = e l6 U (2.30) 
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where ca G C are functions of the coupling constants of the potential V, and the sum goes 
over all irreducible polynomial representations of U (N) which may be parameterized by 
their highest weight components n\ > . . . > un > associated with the lengths of the 
lines in the corresponding Young tableaux. Since 



Xn(U) 



Xn(U 



where 



Ci(n) 



N 

Y, n k 

k=l 



(2.31) 



(2.32) 



is the linear Casimir of the U(N) representation n (the total number of boxes in the Young 
tableau), it follows that the constraint on the SU(N) integration in Q2.29|) is a restriction 
to terms with linear U(N) Casimir C\ equal to iV 2 . We may therefore write 



k N 
N 



[dU ] e 



NtrV(U ) 



SU(N) 



(2.33) 



Ci=JV 2 



where the bracket means to restrict the character expansion of e NtrV(Uo) to those n with 
Ci(n) = N 2 . The role of the logarithmic interaction in (2^28 ) may therefore be charac- 
terized as restricting the matrix integral to completely filled Young diagrams. 

The SU (N) integral in fl2.33p can be evaluated using standard techniques from lattice 
gauge theory [pOfl. For this, we introduce the generating function 



W(J) 



[dU ] e 



tr JU 



SU(N) 



with the property 



[dU ] f(U ) = f (Jj) W(J) 



SU(N) 



J=0 



(2.34) 



(2.35) 



where J is an arbitrary N x N matrix and / an arbitrary function on SU(N). Using the 
invariance of the Haar measure [dU ] under SU(N) rotations, we have W(J) = W(UJV), 
U, V G SU(N), from which it can be shown that the generating function may be expanded 
in powers of determinants of J as 



W(J) 



(27T 



-1 oo 



N 



E 



(det J) 1 



Applying this result to the partition function ( p. 33 ), the restriction C\(n) 
that only the k = N term of the expansion (|2.36|) contributes, giving 



(2.36) 



iV 2 implies 



j=o 



(2.37) 



14 



Comparing (|2.3| ) and ( 2.37) , we see that the finite N expressions for the fermionic partition 
function in the Hermitian two-matrix and unitary one-matrix representations are dual to 
each other. We can make this duality precise by writing (|2.37|) as 



(27T) 



iV 2 



dJ 



^> det^J 



gl(N) 



dJ 



tet N (-i e 



i tr <f> J 



gi(N) 



d4>8{(j>)de% N {-i-^) e NtrV ^ = det 7V (^) e NtlV ^ 



gl(N) 



<t>=o 



(2.38) 



Some further aspects of the unitary representation of the fermionic matrix model are 
discussed in appendix A. In section 5 we will show explicitly that the equations of motion 
of the unitary and fermionic one-matrix models are identical. 



3 Orthogonal Polynomial Solution 

In this section we shall demonstrate how to generalize the orthogonal polynomial technique 
[ pT| to the adjoint fermion one- matrix model ( |1 . 1| ) , with the goal of solving the matrix 
model in the large N limit. Although the fermionic matrix model is not naturally an 
eigenvalue model, since it is not possible to diagonalize a Grassmann valued matrix, 
we can exploit its equivalence with the unitary one- matrix model ( |2.28|) . We will thereby 
define the orthogonal polynomials of the fermionic matrix model ( |1.1| ) to be the orthogonal 
polynomials associated with the corresponding unitary eigenvalue model. This definition 
will lead, as we shall see, to a well-defined real-valued solution for the fermionic matrix 
model. 



3.1 General Properties 



Let us consider the fermionic partition function in the representation ( |2.20|J2.2T| ), which 
involves the measure 



M z ) = ~^n+t dz ( 31 ) 

on the unit circle zz — 1. The partition function ( [2.20|) may then be written as the 
determinant of the moment matrix of the measure ( |3 . 1| ) , 



djq det 



(3.2) 



where d^ is an irrelevant numerical constant and we have defined the inner product 



F(z) 



G(z)) = I dfi{z) F(z)G(z) 

2=0 



(3.3) 
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on the space C[z] <S> C[z] of Laurent polynomials on the unit circle. The construction of 
orthogonal polynomials to evaluate inner products on the circle such as ( |3.3|) with real, 



positive definite measure dfx was discussed long ago in the mathematics literature [22 



and subsequently in the physics literature |23| . Some aspects of generic unitary matrix 
models are also dealt with in |24|, p5| . In the present case, the measure (|3.1|) is complex- 
valued, and we must deal accordingly with defining the appropriate system of orthogonal 
polynomials. 

These are the monic polynomials ($ n (z), A m (z)) of order (n, m) which form a complete 
set in the space C[z] ® C[z] and which are bi-orthogonal in the measure ( |3.1| ). 

$ n (z) | A m (») = h n 5 nm (3.4) 

where h n are some constants. They are normalized as 

n— 1 m— 1 

$ n (z) =Z n -J2 Pn,k Z h , K m {z) =Z- m -J2 W Z~ k (3.5) 



fc=0 k=0 



where the coefficients p n) k and can be formally obtained from the usual Gram-Schmidt 
orthogonalization procedure applied to the monomials (z n , z~ m ) and the inner product 
(|3.3|) . In fact, by iterating ( |3.5|) we find that the monomials (z n ,z~ m ) can be expanded 

as 





l n-1 fei-l 


fcj-_l-l 


Z* 1 = *»(*)+ E| 


E E-- 


■ E PnMPkiM 


3=0 ' 


yfci=0 fc 2 =0 


kj=0 


m 


/ m— 1 fei— 1 


fej_i-l 


z~ m = A m (^) + E 


E E- 




3=0 


\fc 1= fc 2 =0 


k j= 



h,, ,/, A, f:)J (3.6) 

It follows that any polynomial in the space C[z] ® C[z] of degree (n, m) can be expressed 
as a linear combination of ($^(z), A^z)) with k < n and £ < m. From fl3.5| ) it follows they 
define a change of basis (z fe-1 , z~ £+1 ) i— > (§k-i( z ), ^e-i( z )) of the vector space C[z] ® C[z] 
that leaves the partition function ( |3.2|) invariant and at the same time diagonalizes the 
inner product ( |3.3| ). In particular, from ( |3.4| ) we have 



(In det 



iV-l 



$<_i(z) Aj-i(z)) =d N h%l[R%- n (3.7) 



n=l 



where we have introduced the recursion coefficients 



Rn = P~ (3.* 



and h = § z=0 dfi(z) is the normalization of the measure Q3.1| ). 

The reason why two independent sets of functions corresponding to clockwise and anti- 
clockwise rotations on the circle are required here is the complexity of the measure ( |3 . 1| ) . 
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In the classical case whereby the measure dp, is real and positive |22], the polynomials 
$> n (z) and A n (z) are related to each other by complex conjugation. The present system 
of polynomials may be thought of as a deformation of those for positive definite measures 
on the unit circle. This deformation is essentially encoded in the coefficients p n +i,o and 
/ n +i,o- I n order for the Gram-Schmidt process to work for the complex measure (|3.1|) , 
one needs polynomials in both z and 1/z and a priori these polynomials are unrelated. 
Alternatively, we can think of this doubling as the usual doubling of degrees of freedom 
due to the fermionic nature of the original matrix model. As the orthogonal polynomials 
define a complete bi-orthogonal system of functions, they satisfy the completeness relation 

g $ n (z) A n (z) = = zN+l e _ NV(z) §{z) (3 g) 

71=0 ^ n 

on the vector space C[z] ® C[z], where the Dirac delta-function corresponding to the 
measure dp is defined so that 

dp{z) F{z) 5^{z - z') = F{z') (3.10) 

2=0 

for any function F(z) on the unit circle. The relation (|3.9|) is understood as an analytical 
continuation (to the distribution space that is the completion of C[z] <8> C[z]) which can 
be thought of as a representation of the potential Vfyip) in terms of the orthogonal 
polynomials. 

The problem of evaluating the partition function of the fermionic matrix model there- 
fore boils down to evaluating the coefficients R n appearing in (|3.7[). For this, we will 
first derive a few properties of the system of orthogonal polynomials introduced above. 
Completeness implies the relations 

n -i m 

z$ n (z)=$ n+1 (z) + J2'P { k n) ®k(z) , -A m ( 2 )=A ra+1 ( Z ) + ^4 m) A t ( Z ) (3.11) 

The coefficients V^jf 1 and C™^ completely determine the form of the orthogonal polyno- 
mials and hence the partition function. To see this, we use ( |3.5| ) and ( |3.11 ) to write 

n-l 

d> n (z) = z^ 1 {z)~Y,^k~ 1) ^k{z) 



k=0 

n-l n—2 k—l 

k=l k=0 j=0 



Comparing the various polynomial coefficients in ( |3.12| ) with those of the definition ( |3.5| ) , 
we arrive at the iterative relations 

Pn,n-l = Pn-l,n~2 + (3-13) 

n-l 

Pn,k = +Pn-i,k-i+ E P5,hV { r 1] , l<k<n-2 (3.14) 



j=k+i 



n-l 

Pn,o = V^ + Em^ (3-15) 
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The relations (|3.13|) and (|3.15|) can be iterated straightforwardly and by induction we 
have 



n-l 
j=0 



(3.16) 



Pnfl 



n-3 / n-l Jfc-i^ 1 

S n) + E E E ••■ E L) ---?f ; '^r 

&=i \ii=i J2=i ife=i 



-1 il-l jn-3-l 



31=1 J2 = l jn-2 = l 



(3.17) 



which can be substituted into (|3.14j) to iteratively determine the remaining p n> k- Analo- 



gous relations exist between the coefficients and C k . 



We will also need higher degree versions of the completeness relations ( p. 11 ). For this, 
we define constants P^\ by 



n+l 



z e <S> n (z) = £ $ k (z) with PH +e = 1 



[(} 



(3.18) 



k=0 



and use (|3.11|) to write 



n+t+l n+l I k 



fc=0 



fc=0 



3=0 



Using completeness of the orthogonal polynomials we arrive at a recursive relation for the 



coefficients P^ k , 



n+l 

p[*+-i] = pM , v P M -P (j) 

n,k n,k—l ' / , n.j ' k 



<k <n + £+l 



j=k 



with PjfLi = 0, P^™'' = for > n, and the initial condition P^ k = 5 nfc . 



>(«) 



(3.20) 



3.2 Recursion Relations 

We will now determine the coefficients appearing in the completeness relations using 
the above properties of the orthogonal polynomials. This will lead to a set of recursion 
equations for the coefficients R n which thereby completely determines the solution of 
adjoint fermion one-matrix model. Consider the bi-orthogonality relation 



$n+l(2) 







< m < n 



(3.21) 



Using the completeness relations ( |3.11| ) we can write ( |3.21| ) as 



z® n {z)\k m {z))=V ( £ ) h r 



(3.22) 
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The inner product in G3.22j ) can be evaluated by grouping the factor of z with the poly- 
nomial A m (z) with the result 

| A m (z)) = {$ n (z) | A m ^(z) + 0(z- m+2 )-l mfi z) = -l m>0 (z^ n (z) | l) (3.23) 

where we have used the fact that, aside from the term / m z which comes from the constant 
part of the A polynomials in (|3.5|) , the function z A m (z) is a sum of polynomials of degree 
< m — 1. The final inner product in ( |3.23| ) can be evaluated in a straightforward fashion 
using the relation 



z&Jz) 



A n +i(z 



h 



n+l 



(3.24) 



which follows from ( |3.11|) . Again, by grouping the factor of z with A n+ i(z) one finds 



z$Jz) 



h„ — h 



n+l 



I 



n+1,0 



(3.25) 



Substituting ( |3.25| ) and ( |3.23| ) into ( |3.22| ), we find a concise relation between the constant 
terms of the $ polynomials and the completeness coefficients Vx? , 



lm,0 h n+ i — h T 



I 



n+1,0 



kr. 



Q m [R n+ i-Vj Yl QkRk for 0<m<n 

k=m+l 



(3.26) 



Q n (R 



n+l 



where we have defined 



r 



Q, 



for m = n 

ln,0 



I 



n+1,0 



(3.27) 



The completeness relation (|3.11|) for the $ polynomials can therefore be written as 



z$ n (z) = $ n+ i(z) + 



h n +i ~ h n IkQ 

—, Lrw) 

'n+1,0 fc = o IL k 



(3.28) 



Subtracting this relation from itself under the shift of index n — > n — 1 leads to the 
three-term recursion relation for the $ polynomials 



$„+i(z) = z$ n (z) - Q, 



Rn+l 



Rr, 



zR n $ n Jz 



(3.29) 



Analogous relations for the A polynomials may also be derived. Equating the constant 
terms on both sides of ( |3.29|) yields 



Pn+1,0 _ n Rn+i — 1 



Pnfl 



Rr, 



(3.30) 



so that 



Pm,0 h n+ i — h T 



Pn+1,0 



Rk Rk 



-Qn-rra \ Rrn Rn+l 1 j-j - 1 We - 1 We 

-l - 1 fc=m+1 Qfc ^Rfe+i - l 



1 — Rr, 
Qn 



Qm Rn 

for m = n 



for < m < n 



(3.31) 
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Substituting (|3.31|) into the completeness relation (|3.11|) for the A polynomials and sub- 
tracting the resulting expression from itself under the index shift n — > n — 1, we arrive at 
the three-term recursion relation 



A n (z) = -R n A n _i{z) -Q n [A n+X {z) A n (z 



(3.32) 



The recursion relations (|3.29|) and (|3.32|) have been derived without any reference to 
the particular details of the model, i.e. they hold independently of the precise form of 
the potential V of the fermionic matrix model. Consequently, these relations are merely 
'kinematical' in origin. In order to solve for the dynamics of the matrix model we need 
more information. From a more pragmatic point of view, in order to solve for the partition 
function we need to find a solution for the coefficients R n , but the recursion relations also 
involve the quantities Q n . More information is needed to link these two objects. As 
we will see, one only needs two 'dynamical' relations to obtain a closed set of recursion 
relations for any polynomial potential. The first one is given from ( |3.4[ ) and ( |3.5| ) as 

($ n (z) \K l _ l {z)) = {l-n)h n (3.33) 

Integrating by parts on the left-hand side of ( |3.33|) gives 

(N + l)($ n (z) | \ A n _i(z)) - N(V(z) $ n (z) | An-i(z)) - (& n {z) | K-i(z)) = (1 - n)h n 

(3.34) 

Using the completeness relation ( |3.11|) for the A polynomials and the fact that & n (z) = 
n& n -i(z) + 0(z n ~ 2 ), we arrive at 

N(y\z) <f> n (z) | A n -!(z)) = (N + n)h n - nhn-i (3.35) 

For a polynomial potential ( |1.2[ ) of degree K, we may use ( |3.18| ) to write ( |3.35| ) as an 
equation for the recursion coefficients R n , 



Rr, 



n 



+ 



N 



K-l 



n + N n + N 



V" n P [k] 
/ t yfc+1 r n,n-l 



(3.36) 



k=l 



Closely related to the first one, the second 'dynamical' relation is given by the bi- 
orthogonality relation 



o 



Integrating this by parts gives 

N(v'(z)$ n (z) | A n - 2 (z)) 



A n - 2 {z) 



(3.37) 



(3.38) 



The right-hand side of 



can be evaluated using 



and orthogonality, and 



keeping only those terms in the $ polynomials which have a non-vanishing overlap with 
A n _ 2 {z) in (|3T38p . This gives 



n z 



n—l 



— (n — \)p. 



n,n—l 



+ 0(z n - 3 ) \A n _ 2 {z)) 



[n - l)p n ,n-i ~ np 



n-l,n-2 



h 



n-2 



(3.39) 
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By defining 



Pn,n— 1 

N 



we thus find that the second 'dynamical' relation reads 



nP n -i - (n - l)P n = 9k+i Pll 



K-l 

E 

k=l 



>[*•■] 



,,n-2 



(3.40) 



(3.41) 



The coefficients P n can be related to the coefficients Q n and i? n by using (|3.16f) and fl3.26|) 
to get 

N(P n+1 - P n ) = Q n (R n+1 - 1) (3.42) 



To summarize the results of this subsection, we note that the right-hand sides of ( p.36| ) 
and (|3.41| ) can be evaluated, using the recursion formula ( p.20| ), in terms of the complete- 
ness coefficients Vj$. Using ( |3.26|) , it follows that, for any polynomial potential, (|3.36| ), 
(|3.41| ) and ( 3.42|) form a closed set of recursion relations involving only the coefficients 
R n , Q n and P n . An explicit system of fermionic orthogonal polynomials is constructed in 
appendix B. 



3.3 Free Energy 



The final quantity we need for the solution of the fermionic matrix model is the normal- 
ization ho of the measure (|3.1|) . This integral can be evaluated explicitly to give 



h 



dz 



N 



^NV(z) 



2 = 



z N+l N \ Q Z N 



2 = 



(3.43) 



We substitute these quantities into ( |3.7| ) and normalize by the Gaussian partition func- 
tion Z^ uss = (-l)Ma (N 9i ) n2 for which V{z) = g t z, R n = n/{n + N) and h = 
2m\Ngi) N /N\. For a polynomial potential, we may evaluate (|3.43|) explicitly using ( 2.111) , 
and using the multinomial theorem we thus find that the free energy Fn = 4? log „^ ss 
is given by 



'JV 



E(-i) n (^-i)' E 



E 



n 



n=l 

K [N]j 

x n n 

j=2 e j= i 

1 N-l 



0<fc f 



(2) ,„(2) 



■.« [JV ] 2 ^ n 



0<fc, 



He,m k 



K 

n 

m=2 



N 



n=l 



n\, (n + N)R r: 

— log 

N J n 



(3.44) 



where the coefficients R n are determined from (FT 



Generally, in the large- N limit, the variable x = n/N becomes a continuous parameter 
in the interval [0, 1]. We assume that in this limit the recursion coefficients TZ n (any of the 



21 



quantities R n , Q n or P n ), determined as the solutions of the recursion equations derived 
in the previous subsection in the large-iV limit, become continuous functions TZ(x) of 
x G [0,1]. This is justified by the form of the recursion relations which suggest the 
replacement, TZ n — ► = TZ(x). It follows that shifts in the index n are related to 

sub- leading contributions in the large N limit, 

R„ + ^ + - I - (-) (3.45) 



and the second, finite sum in fl3.44|) can be approximated by a simple one-dimensional 



integral over x G [0, 1] at N — oo. To deal with the first term in ( p.44 ), we need to 



determine the large N limit of the normalization constant h . For this, we note that there 
is nothing particularly special about the choice of integration contour in ( |3.43| ), as it was 



merely introduced in (|2.19| ) as a means of evaluating the derivatives in (|2.9| ). In particular, 
since its integrand is an analytic function, we can deform the contour arbitrarily in the 
complex plane, and in particular to one out at infinity. In the large N limit, the one- 
dimensional integral ho may then be evaluated using the saddle-point approximation. In 
this way we find that the genus free energy is given by 

F = lim F N = f + fdxCL- x) log (l±l¥M ( 3 . 46 ) 



where 

/ = nC)-log<?iC (3-47) 
with ( the solution of the order K algebraic saddle-point equation 

CnO = l (3-48) 

for the branch which has the perturbative Gaussian limit Clg2=...=3if=o = l/<7i and for 
which the function ( |3.47| ) is real-valued. 



4 Topological Expansion 

In this section we will solve for the continuum limit of the discretized random surface 
theory represented by the adjoint fermion one-matrix model. We will consider two limits 
of the matrix model. The first one is the naive large- iV limit which captures the lead- 
ing critical behaviour and gives a non-perturbative solution to the problem of counting 
fermionic ribbon graphs of spherical topology. The second one is the double-scaling limit 
in which an appropriate coupling constant is tuned in the limit iV —>■ oo in order to obtain 
contributions from all orders of the ■h expansion of the free energy and which yields the 
complete topological expansion of the fermionic matrix model. 
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4.1 Critical Behaviour 



We will start by making some general remarks concerning the large N limit of the fermionic 
matrix model. From ( 3.46 ) we see that there are two contributions to the N = oo free 
energy. These two quantities represent very different phases of the system. The first one 
is given by ( |3.47| ) which is determined by the saddle-point equation ( |3.48| ). This equation 
is immediately recognized as the critical equation describing the continuum limit of an 
abstract branched polymer || [F]. In fact, the explicit derivative expression for the measure 
normalization ( [3.43 ) coincides exactly with the partition function of the N dimensional 
fermionic vector model with the same polynomial potential V ||14|| . The latter quantity is 
related to the generating function for branched polymer networks by a simple mapping 
of its Feynman diagrams and an analytical continuation N i— ► — ~ of the fermionic vector 
dimension. The analytic reason for the appearance of branched polymer behaviour in the 
fermionic matrix model is that the normalization of its measure is not of sub-leading order 
in N. From a graphical point of view, its appearance is clear from the analysis of section 
2.2. The entropy factors associated with the number of Feynman graphs is drastically 
reduced in the fermionic case because of the cancellations which occur between twisted 
diagrams. There is a class of diagrams with no twists in the fermionic matrix model, and 
these have a tree- like growth in the large N limit. The reduced entropy of the remaining 
twisted graphs is then comparable to that of the graphs which produce a polymer-like 
behaviour. This feature is unique to fermionic matrices, and it is the property which 
enables the construction of a model of branched polymers using supersymmetric matrix 



models [13] via the coupling of the fermionic matrix model to an ordinary, complex matrix 
model which has the effect of cancelling the set of twisted diagrams. In fact, it is precisely 
the twisting mechanism described in section 2.2 that enables one to isolate graphs with 
tree-like growth in the fermionic case. The untwisted diagrams may then be mapped 
onto branched polymers similarly to the case of the cactus diagrams which appear in 
supersymmetric matrix models |13| . 

The twisted graphs are associated with the second contribution in fl3.46| ) and, as we 
will demonstrate, they lead to the usual surface effects in the theory. Thus the orthogo- 
nal polynomial formalism that we have developed naturally splits the generating function 
into a piece corresponding to the tree-like graphs and a piece corresponding to the twisted 
diagrams. The continuum limit of the latter ensemble of graphs describes a surface theory 
and describes an effect which is of order iV 2 , while the former ensemble which has a poly- 
mer growth produces an effect of order N. One typical feature of the polymer generating 
function is that it becomes non- analytic in the large N limit at the Gaussian point gt = 0, 
k > 2. However, this effect is of sub-leading order in N, so that a large N analysis with 
arbitrary coupling constants will miss the branched polymer phase transition. There is 
a "barrier" at the Gaussian point which separates the surface-like continuum limit from 
the branched polymer behaviour. From a geometrical point of view, the fermionic ma- 
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trix model probes an intermediate lattice phase where the number of Feynman diagrams 
contributing to a given process is greater than that expected of a vector model but less 
than that of a matrix model. Thus the matrix integral generates a random surface theory 
which contains an "internal" branched polymer phase. 

This decrease in the number of graphs as compared to the usual matrix models will 
be responsible for the Borel summability of the all genus expansion of the free energy of 
the fermionic matrix model that we will prove in this section. As discussed above, the 
essential critical behaviour of the fermionic matrix model is encoded in the large- N limit 
of the recursion coefficients R n . If we introduce the quantities r n defined by 

K = ^ (4.1) 

with r = cIn as in (|3.2|) , then the recursion coefficients are given by R n = r n+ i r n _i/r^ 
and the free energy by Fn = Fn, where T n = log(r n /r^ auss ). Replacing sequences by 
functions of x e [0, 1] as prescribed by ( |3.45| ), we then have to leading order in N, 

X ox 2 

The relations ( |4.2j ) show that the function R(x) is related to the specific heat uq of 
the matrix model, and also that any singularity in the free energy will occur at the 
boundary x — 1 of the domain of R(x). These facts will enable us to construct the 
topological expansion of the fermionic matrix model from the recursion coefficients i? n 
of the orthogonal polynomials. We remark that it is possible to demonstrate from the 
constrained sum over SU(N) group characters of section 2.3 that the partition function 
is real and has a definite value for any N. Combining this with our knowledge from large 
iV loop equation calculations, wherein the system goes to a definite limit, one can argue 
that there is a definitive expansion of the free energy. The phase transition at large 
iV can thereby be understood as a sort of percolation transition in the Young tableaux, 
whereby the system evolves into a phase in which the number of filled SU(N) Young 
tableaux grows asymptotically. 



4.2 Main Recursion Relation 



We will now begin quantifying the above discussion. For definiteness, in the remainder of 
this section we will deal primarily with the simplest non-Gaussian model which is given 
by the quadratic potential 

V(z)=z+ 9 -z* (4.3) 

although, as we will discuss, the results we obtain are universal. The potential ( fO| ) is 
difficult to treat within the loop equation approach because the endpoints of the support 
of the corresponding spectral density are located asymmetrically in the complex plane 
[fTD| . Indeed, one power of the technique that we develop is that it is insensitive to the 
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parity symmetry (or lack thereof) of the potential. This will account for the universality 
properties of the solution that we shall find. In this subsection we will derive the main 
recursion relation corresponding to the potential ( |4.3|) that will be used to construct the 
topological expansion. 

To solve the recursion relations ( |3.36| ), ( |3.41| ) and ( |3.42| ) corresponding to 
first use ( |3.2(J| ) and ( |3.26| ) to determine the coefficients 

Qn-lQnRn{Rn+l ~ 1) 



p[l] 
r n,n— 1 



P, 



[1] 



71,71 — 2 



71-1 

<r>(«) 
'ti-2 



Qn~2Qn~lQnRn-lRn(Rn+l ~ 1, 



This leads to the closed set of recursion relations 

(N + n)Rn-n = g N Q n ^Q n R n (R n+1 - 1) 

nP n -l - (n - l)P n = gQn-2Qn-lQnRn-lRn{Rn+l-^) 

N(P n+1 -P n ) = Q n (R n+1 -l) 
From the relations ( |4.6| ) and ( |4.7p one can solve for the coefficients P n as 

NP n = g N Q n _ 2 Qn-lQnRn-lRn{Rn+l - 1) +TlQ n „i(R n - 1) 



^|), we 

(4.4) 

(4.5) 
(4.6) 
(4.7) 

(4.8) 



Substituting ( fOf ) into ( |4.7| ) and using ( |4.5| ) to reduce the degree of the Q n 's in the 
resulting equation yields 







n 



Qn-tRn [(N + n + l)R n+1 - n - lj - Q n - 2 Rn-\ [(N + n)R r , 
+ nQ n (R n+1 - 1) -nQ n -i(Rn - 1) 

We now use ( f4.5| ) to define the quantity 

Qim - nn _ 1 (N + n)R n -n 
n n {R} = QnQn-l ~ ~ NRn{Rn+i _ X) 

and multiply the relation fl4.9|) through by Q n . Using ([4.6|) along with ( |4.8| ), we can then 
solve for Q n in terms of the R n 's alone as 



(4.9) 



(4.10) 



Q n {R} 2 = n n {Ry 



n(R n - 1) - Rn 


(N + n + l)R n+ i - n - 1 




nn n {R}(R n+1 -l)- 


&<n-l{R} Rn-l 


(N + n)R n - n 



(4.11) 



We may therefore write down a recursion relation involving only the coefficients R n in 
the form 

Qn{R} 2 QnMR} 2 = ^n{R} 2 (4.12) 

where the function Q n {R} 2 is given by and tt n {R} by (glQ|) . Notice that the 

coupling constant g completely cancels out in this final equation which is an involved non- 
linear recursion relation for the i? n 's solely in terms of n and N. The dependence of R n 
on g comes about purely as a boundary condition and can only be recovered after solving 
for the dynamics of R n . However, we will see that the entire topological expansion can 
be constructed without ever knowing this dependence explicitly. All information about 
the continuum limit will come from the details of the main recursion relation ( [4.1 2| ) . 
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4.3 Planar Limit 



The spherical continuum limit of the matrix model is found by taking the limit N — ► oo. 
Let us first describe the polymer contribution f Q in (|3.46|) . For the potential (fOj ), the 
solution of the quadratic saddle-point equation ( |3.48|) which is regular at g = is given 
by 

C=^(v / l+^-l) (4-13) 
so that the corresponding fermionic vector free energy (|3.47 ) is 

1 1 / / \ , 1 



/o(^) = 2 + ^ (yi + ^-lj-log- (yi + 40-lj (4.14) 

There is a phase transition at g — g^ = — | below which the free energy becomes complex- 
valued. Near this critical point, (|4.14 ) behaves as fo(g) ~ (g — g^) 3 ^ 2 , which identifies 
the corresponding string susceptibility exponent as 7^ = +§. The phase transition is 
therefore associated with the evolution of the system into a pure branched polymer phase 
of two-dimensional quantum gravity. In fact, the function fl4.14j) is related to the usual 
bosonic vector model free energy T veci which is the generating function for pure, connected 



branched polymer graphs, by fo = 1 — 2jF vec WM. The iV dimensional Euclidean radial 
coordinate r of the bosonic vector model is related to the saddle-point (|4.13|) by r 2 = | (. 
One can also compute the one-loop fluctuations around the saddle-point (|4.13|) and show 



that the resulting free energy is logarithmically divergent at g = JT4j. Therefore, in 
the regime g > there is no polymer behaviour and one may imagine a "barrier" in the 
large N limit of the theory at the Gaussian point g = 0. It is further possible to carry 
out a double scaling expansion of the fermionic vector model partition function (|3.43|) in 
the limit iV — > 00, g — > g^ with the parameter N(g — g(°)) 3 / 2 held fixed. The resulting 
expression is an alternating, Borel summable series which can be expressed as a Bessel 
function @. 

We now consider the large N behaviour associated with the second term in the free 
energy ( |3.46|) . Replacing the discrete R n coefficients with the large- N continuous function 
R(x) as prescribed by ( |3.45|) , after some algebra we obtain from the recursion relation 
( |4 . 1 2| ) , at leading order in N, the first order non-linear ordinary differential equation 

dR r(R-1) 2 (x + 2xR-3(1 + x)R 2 ) 

dx ~ x 2 -2x 2 R + Ax{l+x)R 2 -2(l+x){2 + 3x)R 3 + 3{l+x) 2 R 4 ^' ' 

The dependence of R on the coupling constant g arises only from the constant of inte- 
gration of (|4.15 ). It is known from the general theory of ordinary differential equations 
26] that the solutions R(x) to equations such as ( |4.15|) possess algebraic non-analytic 



behaviour. Since dR/dx in ( [4.15| ) is given as the quotient of a quintic polynomial in R 
by a quartic one, it follows that the function R(x) is finite for all finite values of x [p6 



The only singularities which can arise from this differential equation occur when the de- 
nominator on the right-hand side vanishes and the numerator is non-zero. In light of the 
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discussion of the previous subsection, we will demand that these singular points occur at 
the boundary x = 1 of the domain of the function R(x). At x = 1, the denominator of 
the right-hand side of ( |4.15 ) has four distinct zeros R c determined as the solutions of the 



quartic equation 

1 - 2R C + 8R 2 C - 20R S C + Y1R\ = (4.16) 

Denoting the corresponding critical value of the coupling constant g by g c , we make 
an ansatz for the form of the function R(x; g) near the critical point, 

R{x- g) = R c + a-(g c - gx)~^ + ... (4.17) 

where a is some constant and the ellipsis denotes terms which are less singular at the 
critical point. This ansatz ensures that for g ~ g c , the singularities occur at z = 1, 
or alternatively that at x = 1 the free energy T{x) becomes non-analytic at g = g c . 
Substituting ( [4.17|) into the differential equation ft4.15|) , we find for each branch R c of the 
equation ( |4.16| ) the leading behaviour 

dR 1 , A , 

— ~ (4.18) 

dx R c -R{x) v ; 

Comparing with ( |4.17[ ) fixes the exponent 7 str = — |, and so the critical point g = g c of the 



fermionic matrix model describes a continuum limit that lies in the same universality class 
as pure two-dimensional quantum gravity in the planar limit. Note that this criticality 
argument requires no knowledge of the precise value of the critical coupling g c . The same 
will be true of the double scaling limit which will be analysed in the next subsection. It 
is possible to show using loop equations |§ that g c > 0. Thus the continuum surface 
behaviour occurs on the opposite side of the Gaussian point relative to the branched 
polymer phase transition, and the two phases cannot be connected together in a smooth 
way. 

We conclude our analysis of the planar limit of the adjoint fermion one-matrix model 
by briefly describing how it extends to more general potentials. For example, consider a 
potential of the generic form V(z) = z + j^z K , K > 2. The recursion equations ( |3.36| ) 



and ( p.41| ) in this case require knowledge of the coefficients P^-i an d determined 
from 2 Up and ( |3.26| ). It is straightforward to see from these relations that P^-i ls °f 
degree K in both the QnS and the i? n 's, while -P^L 2 is of degree K + 1. In the large N 
limit, the recursion equations will therefore assume the form 

(l + x)R-x = gQ K Wi[R] (4.19) 
BP 

P- x °— = gQ K ^W 2 [R] (4.20) 

ap 

— = Q(R-l) (4.21) 

where and W^t-R] are polynomials of degree K and K + 1 in the function R, 

respectively, which each contain the factor R — 1. Using ( |4.19| ) we can solve for Q K and 
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write ( f4.20|) as an equation determining the function Q. Comparing with the solution for Q 
obtained using ( 4.21|) , the relation ( 4.2C ) then becomes a first order, linear, homogeneous 
differential equation for the function P. Substituting the resulting solution for P back 
into ( |4.19| ) yields an integral equation for the function R which, upon differentiation, can 
be transformed into a first order non-linear ordinary differential equation of the form 

dR = J 2 [R; x] ^ 22 ^ 

dx J\ [R; x] 

where Ji[R;x] and Ja[i?;x] are polynomial functions of their arguments of respective 
degrees 2K and 2K + 1 in R which are independent of the coupling constant g. Arguing 
as we did above, the differential equation ( |4.22[ ) will admit critical behaviour for the 
function R of the square root type, i.e. that of pure gravity. Provided that the zeros of 
the function Ji[R;x] are non-degenerate at x — 1, this will be the only type of critical 
behaviour that the model admits. It is of course natural that pure gravity exists for a 
generic matrix potential. Multicritical points require more complicated potentials with 
two or more coupling constants and an appropriate fine-tuning of them. In this case the 
differential equations derived analogously to those above become highly non-linear and 
quite involved. But the general picture will be the same, namely the function R will be 
determined by some differential equation whose polynomial coefficients (in R and x) can 
be tuned to obtain higher order critical points. In this way we can recover, in the planar 
limit, the standard universality classes of conformal matter coupled to two-dimensional 
quantum gravity [£§]. 



4.4 Double Scaling Limit 

In this subsection we will study the all genus expansion of the specific heat 

oo 

u(g, N) = J2 N- 2h u h (g) , u h (g) = F'^g) (4.23) 

h=0 

where F h (g) is the genus h contribution to the free energy F N (g) and the genus h suscepti- 
bility is a homogeneous function of degree — 2h, Uh{a- g) = a~ 2h - Uh{g)- From the analysis 
of the previous subsection, we know that the leading singular behaviour of ( |4.23| ) at genus 
zero is Uo(g) ~ \Jg — g c . In order to keep the partition function Z (g) ~ e N2p o(9) finite in 
the large N limit, we should therefore perform the double scaling limit g — > g c , N — ► oo 
by keeping fixed the parameter 

n = N^{g c -g) (4.24) 

We now define a scaling function u(k) which captures the contributions to the string 
susceptibility from all genera in the double scaling limit g — > g c , N — ► oo with the variable 
( ^4.24| ) fixed. For this, we introduce the large- N expansion parameter e = and write the 
total susceptibility (|4.23|) in the vicinity of the critical point as 

u(g, N) = (e 2 ) 4 / 5 a + e 2 u{k) + (e 2 ) 6 / 5 a X K + ((e 2 ) 7 / 5 , k 2 ) (4.25) 
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where constants. By defining the new effective variable 

t = e-^{g c -gx) (4.26) 

in the large N limit, this scaling function may be determined from the £ = k limit of a 
scaling ansatz for the function R(x) in terms of an unknown function 

R(x) = R C + e 2/5 u(£) (4.27) 



where R c is a solution of (|4.16|) . By substituting (|4.27|) into fl3.46| ) and changing variables 



from x to £ in the integral, we find that, up to irrelevant constants (which may be absorbed 
by suitable rescaling using the homogeneity of the specific heat) and terms which vanish 
as e — > 0, the free energy in the double scaling limit is determined as 



F(«)=/ + / d£(K-0u(0 (4.28) 



-4/5 



} 9c 

It follows that 

u{k) = F"(k) (4.29) 

and so the problem of obtaining the topological expansion of the fermionic matrix model 
is thereby reduced to the task of finding the solution of the main recursion relation ( |4.12[ ) 



for the function R(x) in this special limit. Note that the polymer free energy ( 4.14 ) 
contributes only an irrelevant constant in the double scaling limit about the positive- 
valued critical point g c . 



We now rewrite ( |4.12j ) using ( |3.45| ) and take the limit of large- N while holding the 



quantity k fixed through the relation g = g c — kN 4//5 . We then substitute in the scaling 
ansatz ( 4.27 ) and rewrite derivatives according to the change of variables (|4.26|) , i.e. 

d x e- 4 / 5 <%. After some algebra, we find, as the coefficient of the leading order e ' 

term, a non-linear differential equation for the specific heat 

= 6g c (3 - 22R C + 80R 2 C - 192i2jj + 376R 4 C - 632R 5 C + 752R 6 C - h\2R 7 c + 144^) 
x u{i) u\i) + R c (l - 5R C + 10R 2 C - 10R 3 C + ARf) 

x QR C (R C -if (6R 2 C -2R c -l)-gl (l-2R c + 20R 2 C - 56i^ + 36i? 4 ) u M (£) 

(4.30) 

We now integrate ( |4.30|) up once and use an inessential shift of the independent variable 
£ to eliminate the constant term. After applying ( J4.16p and a rescaling using the homo- 
geneity of the function (which preserves the relationship ( [4.29] )), we arrive at the 
parameter-free equation 

u"(0+u(0 2 = Z (4-31) 

The non-linear differential equation (|4.31| ) governs the behaviour of the partition function 
of the fermionic random surface model to all orders in the genus expansion, with the 
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parameter £ 5 / 4 identified as the renormalized string coupling constant. We note again 
the remarkable feature that one never needs to know the precise location of the critical 
point to arrive at this equation. In the present case the critical coupling constant g c can 
be consistently rescaled out of the pertinent equations. This parametric independence is 
indicative of the universality of the double scaling equation ( |4.31|) for the given class of 
generic polynomial potentials. 



The differential equation ( j4.31|) is known as the Painleve I equation and it is solved 



by the first Painleve transcendent. This equation differs from the versions which usually 
appear in matrix models 0] by a change in sign of the specific heat u — > —u. The 
boundary conditions which satisfies follow from the analysis of the planar limit of 
the previous subsection. The genus zero contribution arises in the limit of large, positive 
£ where the leading behaviour is w(£) 2 = £• As a consequence, we are led to postulate an 
asymptotic expansion for tt(£) of the form 

oo \ 
-5fc/2 

•*« S 

k=l 



By substituting ( |4.32j ) into (|4.31|) , the first few coefficients are easily calculated to be 

1 



Ul = 


8 




49 


u 2 = 


128 




1225 


U3 = 


256 




4412401 


W 4 = 


32768 




220680075 


U 5 = 


32768 




2207064977649 


U 6 = 





(4.33) 



and in general they are determined by the recursive equation 



(5fc-6)(5A;-4) 1 k ^ , ^ QA , 

Uk = Uk-i - - V u k u k _ n , k > 2 (4.34) 

8 2 n=l 



With the normalization of the spherical specific heat taken as Uq — 1, we see from fl4.33|) 



that the coefficients of the genus expansion have precisely the same absolute numerical 
values as those obtained for pure two-dimensional quantum gravity, i.e. from the usual 
Painleve expansion [^|, In particular, they have the same high rate asymptotic growth 
in magnitude. However, the most notable feature of the coefficients ( [4.33 ) is the fact 



that they oscillate in sign (In the Hermitian cases, all u k are negative). This raises the 
possibility that the asymptotic series solution ( |4.32j ) may in fact be Borel summable. If 
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so, then the free energy would give an unambiguous definition of the genus expansion of 
pure "fermionic" gravity in powers of N~ 2 ~ £~ 5 / 2 . A numerical solution of the equation 
(|4.31 ) is possible and the result is shown in fig. 1. This confirms numerically that the free 
energy exists as a well-defined function and is real-valued on the positive real £-axis. In 
the next subsection we shall prove that the coefficients Uk of the asymptotic series ( |4.32|) 
alternate in sign to arbitrarily large orders of the genus expansion, and, moreover, that 
there is a well-defined Borel resummation of for £ > 0. 
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Figure 1: Plot of the solution tt(£) versus £ of the Painleve equation with the asymptotic 
boundary condition w(£) 2 — £ f or £ — * +oo. The first double pole singularity (correspond- 
ing to a zero of the partition function) occurs on the negative axis at £ ~ —3.1477. 



4.5 Borel Summability 

In this subsection we will argue analytically, along the lines of [|J, that the asymptotic 
series (|4.32| ) defines a unique function For this, we consider the Borel transform 



BW= Em* (4 - 35) 

where the constant (3 will be self-consistently determined by the condition that this series 
has a finite radius of convergence. A solution of the original Painleve equation ( |4.31| ) is 
then given by 



u(£) = v£ \ l-fdte- t B(fC* /2 )j (4-36) 



The crucial issue now is whether or not the integral in (|4.36| ) actually exists. We will show 



that there is a contour running through the Ret > region of the complex t-plane from 
t = to t = oo along which the integral transform ( ^4.36| ) converges. We will thereby argue 
that the Borel resummation ( f4.36p of the specific heat is well-defined and real-valued. 



The function ( |4.35| ) has singularities in the complex s-plane with branch cuts running 
between them and infinity. If some of these singularities lie on the positive s-axis then 
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the integral definition (|4.36|) is ambiguous. We will now extract the large order behaviour 
of the coefficients Uk and argue that the singularities of the Borel transform B(s) all 
lie on the negative s-axis. From ( |4.35| ) and (|4.36| ) it follows that the coefficients of the 
asymptotic expansion ( [4.32|) are given by 



Uk 



s=0 



(4.37) 



The contour of integration in ( 4.3 7| ) can be extended out to infinity around the cuts of 
the Borel transform in the complex s-plane. In the limit k — > oo, the contour integral in 
(|4.37|) is dominated by contributions along the branch cut of B(s) in the complex s-plane 
which begins at the point s = Sq that is closest to the origin. We then have, for large k, 



Uk 



J dt e ~* J -gL Disc B(t p s)~j^ J dt e-* Disc B {t p s) (4.38) 



soft? (so/s) 1 /^ 

where we have defined the discontinuity of the Borel transform 

Disc5(s) = B+(s) -5_(s) , B ± (s) = B(s±iO) 



(4.39) 



for s a point on its cut. 

Going back to (|4.36|) , we see from ( |4.38|) that we should study the functions w±(£) 
corresponding to the Borel transforms on either side of the dominant cut, 

-5/2> 



u±(0= / dte- t B ± (t^C 

(S ? 5 /2)l//3 



(4.40) 



Since each of the functions ( |4.40| ) solves the Painleve equation ( 4.31 ), it is convenient to 
define their sum and difference 



w d (0 =«+(£) "MO = Disc B(0 



; (0 = I + 



which on using (|4.31 ) are seen to obey the coupled system of differential equations 

<(e)+2« d (o« s (o = o 



(4.41) 

(4.42) 
(4.43) 



These equations can be solved in the WKB approximation. In the planar limit £ — ► oo, 
we self-consistently assume that tt d (0 is of sub-leading order in ( |4.43| ) and thereby find, 
to leading order, the solution w s (£) ~ By substituting this into ( |4.42| ) , to leading 
order we have the first order linear differential equation for the function tt d (£), 



<(£) + 2V£ud(0 = o 



(4.44) 
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which can be solved in terms of Bessel functions as Wd(£) = \/£ ^2 (^5^ £ 5 ^ 4 ) • The solution 
for £ — > 00 is therefore given by 

^^j-Weo.^^^^.j) (4.45) 

We are finally ready to estimate the large order behaviour of the coefficients u^. For 
this, it is convenient to change variables in order to bring the asymptotic expansion for 
i -> 00 to u -> through the definition = £ _5/4 . Then, from (^38|)-( ^4lD and (|445| ), 
we find as — > 00, 

* ~ / -£L ^ Md ~/4^«^^-l) (4.46) 



For each choice of sign the integral (|4.46f) gives, up to constant factors, the result 



u k ~ ( * T (2fc - i) for fc -> 00 (4.47) 



The asymptotic estimate (|4.47 ) contains a large amount of information. Going back to 



the Borel transform ( |4.35| ), the large order behaviour of fixes (3 = 2 and gives 

*(«)~E (-§)"*«* ( 4 - 48 ) 



k 



We see therefore that the Borel transform B(s) is well-defined in an open region of the 
complex s-plane and, moreover, that its first singularity appears on the negative real 
s-axis. Assuming that all of its singularities are so restricted, we can expect that the 
asymptotic series ( 4.32|) can be Borel resummed through the integration in ( 4.36Q to 



define a unique real function tt(£) on the positive real £-axis. The results of a numerical 
analysis (fig. 1) support these arguments. 

From this analysis it follows that the double-scaled partition function of the adjoint 
fermion one-matrix model leads to a completely well-defined genus expansion of the cor- 
responding random surface theory. The coefficients ( [4.47] ) have precisely the same large 
order behaviour T(2k — |) as those of pure two-dimensional quantum gravity. In fact, the 
topological expansion of the fermionic matrix model is identical term by term to that of 
the usual Hermitian matrix models. However, the alternating nature of the asymptotic 
series expansion allows one to define the string susceptibility in an unambiguous way 



through its Borel resummation ( 4.36Q , at least for positive values of the string coupling £. 
In this way we may think of the fermionic matrix model as lending a well-defined version 
of the generating functions provided by Hermitian matrix models. It gives an analytic 
solution to the problem of counting fermionic random triangulations on arbitrary genus 
Riemann surfaces, whose asymptotic expansion is Borel summable but otherwise coincides 
with the usual Painleve expansion of two-dimensional quantum gravity. 
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The usual movable, double pole singularities of the Painleve I equation appear as well 
in the present case, but this time they appear for negative values of the cosmological 
constant £. The existence of such poles is inconsistent with the loop equations of two- 
dimensional quantum gravity |2?J . Pole-free solutions of Painleve equations are provided 
by the triply-truncated Boutroux solutions, but these are complex-valued and do not lead 
to physically acceptable free energy functions. However, in the present situation there is 
a unique, well-defined real-valued specific heat on the positive £-axis! which should have 
a meaningful analytical continuation to £ < 0. It provides an unambiguous definition 
of a whole branch of the random surface theory, and, given the relation between the 
double scaling variable £ and the continuum cosmological constant, the continuum limit 
of the discretized model is indeed well-defined. This suggests that by restricting our 
"fermionic gravity" model to £ > 0, we could define a model of two-dimensional quantum 
gravity with the canonical characteristics order by order in the genus expansion, but whose 
full asymptotic series makes perfect sense and provides a well-defined non-perturbative 
definition of the theory. 

In the next and final section we will argue that this is indeed the case. In the usual Her- 
mitian one-matrix models, the problem with the double scaling limit, in which the original 
matrix integral can only be defined by an analytical continuation, can be traced back to 
an instability in the corresponding double-scaled eigenvalue model. The non-compactness 
of the double-scaled eigenvalue space K and the form of the effective potential for the 
eigenvalues demonstrates that the definition of the critical points of the model is unstable 
to the tunneling of eigenvalues into a different configuration p9| . This leads to instanton 
solutions in the single-well eigenvalue model, and thereby explains the complexity of the 
free energy. However, in the present situation the Grassmann matrix integral is completely 
well-defined at finite N, in contrast to the Hermitian cases, and there is no reason a priori 
to expect that the model becomes unstable in any way at large N. Indeed, with the map- 
ping of the fermionic matrix model onto a unitary matrix model, we see that the matrix 
integral is determined by an eigenvalue space whose topology is that of a circle. Being a 
compact space, there is no asymptotic behaviour in the effective eigenvalue action, and 
thus the compactness could have the effect of eliminating the eigenvalue tunneling prob- 
lem. We see then that the unitary representation of the fermionic matrix model provides, 
at least naively, an eigenvalue description of the random surface model which reflects the 
reasons why the matrix integral is superior to its Hermitian counterparts. We shall tackle 
this problem using an operator theoretic approach which gives the fermionic analog of 
the (generalized) KdV hierarchical structure of the Hermitian one- matrix models [§], [BOf . 
The integrable flows that we shall find are similar to those of the bosonic models, but 
with some very important changes. The most important difference will be the absence of 
a translational symmetry in eigenvalue space. This feature can be interpreted as allow- 



2 We thank A. A. Kapaev for pointing out the relevant mathematical literature |28j on this 
point. 
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ing one to restrict the double-scaled specific heat to positive values of the cosmological 
constant £. It may then follow that the partition function of the fermionic matrix model 
serves as an unambiguous definition of that for two-dimensional quantum gravity. 

We stress that the resulting Borel summability of the fermionic case is very different 
from that of Borel summable bosonic models, such as the conformal field theory of the 
Yang-Lee edge singularity in which the lattice system couples to an imaginary magnetic 
field and is therefore non-unitary Q. The coefficients of the genus expansion of the adjoint 
fermion one-matrix model are all real-valued!. The stability of the fermionic model is 
also quite distinct from the stabilizations provided by supersymmetric and stochastic 
quantization methods 132 1, because these latter models violate the KdV flow structure of 
two-dimensional quantum gravity at a non-perturbative level [33]. In the present case, we 
will see that the appropriate integrable hierarchical structure is a deformation of that for 
gravity, and so the KdV flow structure continues to hold in a certain sense. The precise 
meaning of this deformation comes from the worldsheet interpretation of the fermionic 
one-matrix model. This is the model of fermionic gravity that we alluded to in section 2.1, 
wherein "Penner matter states" are located at the vertices of the corresponding surface 
discretizations. 



5 Operator Formalism 

The discussion at the end of the previous section motivates the development of an operator 
approach to analyse more carefully the properties of the adjoint fermion one-matrix model 
within the orthogonal polynomial formalism. The main purpose is two-fold. First of all, 
it will allow us to relate the partition function of the fermionic matrix model to the r- 



function of an integrable hierarchy pOl 134], |15[. The usefulness of this association is that it 



partitions the solutions of the fermionic matrix model into universality classes which are 
related to well-known integrable hierarchies, and thereby formulates the partition function 
in an invariant way as the solution to a set of differential equations, in particular in the 
double scaling limit. Secondly, given this relationship, we will obtain a clear geometric 
picture of the origin of the alternating, Borel summable genus expansion of the model. 
This will allow us to clarify the description of the string susceptibility as a well-defined, 
non-perturbative generating function for the fermionic gravity model, as alluded to at 
the end of the previous section, and at the same time it will present the appropriate 
generalization of the KdV flow structure that characterizes the gravity model, lending a 
complementary characterization to the worldsheet description. In the following we will 
begin by introducing the formalism and deriving the main equations that will be required. 
We will then identify the integrable hierarchy of which the fermionic partition function is 



3 In [31 1 a similar, but non- unitary, alternating genus expansion was used to formulate a model 
of two-dimensional quantum gravity in terms of integrals over the moduli spaces of punctured 
spheres. 
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a T-function, and derive the Virasoro constraints which must be satisfied by the flows. 



5.1 String Equations 



We begin by introducing multiplication operators Q and Q defined on the space C[z](8>C[z] 
by 

(Q$ m )0z) = z^ m (z) , (QA m )(z) = -A m (z) (5.1) 

2/ 



and differentiation operators P and P by 
(P$ m )(z) =#„,(*) 



(P A m )(z) = z' A'Jz 



(5.2) 



The operators Q and Q are related to each other in a very simple way. From ( |5.1| ) it 
follows that 

(Q $ n (z) | QA m (z)) = | A m (z)) (5.3) 

and so Q is an invertible operator with inverse 



Q 1 = Q' 



where the adjoint of any operator O on 



[z] is defined by 



F(z) 



OG)(z)) = ((O^F)(z) G( 



(5.4) 



(5.5) 



By writing ^(z)]' in two different ways as [(Q$ m )(z)]' and $ m (z) + z<&' m (z), we 
can infer the canonical commutation relation 



P,Q 



Similarly, by considering z [- A m (z)]' we arrive at Q , P 
condition ( |5.4j ) implies the canonical commutator 



(5.6) 



1, which with the unitarity 



F , Q 1 



(5.7) 



Geometrically, the operator P is the canonical conjugate of the operator generating clock- 
wise unit shifts on the circle, while P serves as the canonical conjugate of that which gen- 
erates counterclockwise shifts. The operators introduced above are not all independent, 
but are related to each other through a Schwinger- Dyson equation. This follows from the 
inner product 



(N + l)(<S> n (z) | \A m (z)) 



n(V(z) $ r 



^m(z] 



$n(z) 



(5. 
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which implies the operator identity 

P + P t CT 2 = (JV + 1) Q 1 - N V'(Q) (5.9) 
and ( |5.9|) define the string equations of the adjoint fermion 



The equations (|5.6|). 
one-matrix model. 

Let us now describe some basic properties of the operators introduced above. For this, 
we first rewrite their actions on the same basis $ n of the space of polynomials in a single 
variable on the circle. We have 



:q ±x *»)(*) 
(p *»x*) 



E[Q ±X 

Tt>0 

E[p 

Tt>0 

E [ F 



m>0 



(5.10) 



We denote by + the upper triangular part, including the main diagonal, of the discrete 
operator O when represented in the basis of polynomials in ( 5.10Q . Then O = O — 0+ 
is the lower triangular part of O. 



From the definition ( |5.1|) it follows that the operator Q defines a Jacobi matrix, because 

[Q]nm = for m — n>l (5-11) 
In particular, its pure upper triangular part is given by 

(5.12) 

n>0 n>0 

where E„ i?n are the step operators with matrix elements [E niTO ]^ = 5k n &im- The remaining 
matrix elements of Q are given by [Q] nm = in (|3.26| ), along with the recurrence 
relations derived in section 3. Similarly, the operator Q _1 defines a Jacobi matrix in the 
basis ( |5.10 ) because 



Q 



-i 



for n — m > 1 



Since 



K-l [Q~ 1 ] nn _ 1 = ((Q ®n)(z) I An-l(^)) = | \ K-l(z)) = h n 



it follows that its lower triangular part is given by 

Q- 1 — E ^™ ^",n-l 



(5.13) 



(5.14) 



(5.15) 



n>l 



The remaining matrix elements of Q 1 are given by [Q 1 ] nm = in ( |3.31| ). Further- 
more, from the definition (|5.2j) we see that P = P is purely lower triangular with 



-l = n 



(5.16) 
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and since 



'ra+1 



n,n+l 



P 1 $„)(*) A n+1 (2) = $ n (z) z 2 A' i(*) ) = -(n + lK (5.17) 



it follows that = P+ is purely upper triangular 



n,n+l 



P 

n + 1 

Rn+l 



0, with 



(5.1* 



Note that, in this basis, the (n, n — 1) matrix element of the Schwinger- Dyson equation 
(|5.9|) coincides with (|3.36|) . This follows from the structure of the matrix elements 



h n -l [P T Q 

which implies that P^ Q 



n,n—l 



n,n—l 



M n _ 1 (z)) = -(n-l)h n 



{n - l)R n . 



(5.19) 



5.2 Flow Equations 



We will now derive a system of flow equations for the operators above and the partition 
function (|3.7| ), which will enable us to identify the particular integrable hierarchy at work 
here. We shall describe the evolutions with respect to the discrete set of "time" variables 
tk = Ngk/k of the generic potential (|1.2j). Consider the relation 







Ar 



z)) + (£- k $ n (z)\A n - m (z)) (5.20) 



where 1 < m < n and we have used the fact that "^ n dt ™ (Z> is a polynomial of degree at 



dA n _ m (z) 

most n — m — 1. This relation leads immediately to the evolution equation 



dt k 



§ m (z) 



m=0 



(Q fc - $ ») (z) 



(5.21) 



By taking time derivatives of the equations ( J5.1U ) , we then arrive at the discrete operator 
flow equations 



dQ 

dt k 
OP 
dt k 
dQ- 1 
dt k 
dP j 



dt k 



Q,Q fc 
P , Q! \ 
Q 1 , Qt 
P t , Q A 



(5.22) 
(5.23) 
(5.24) 

(5.25) 



The flow equations (|5.22 ) constitute a discrete KP hierarchy |3"4" |. It is straightforward 
to show that these flows are mutually commutative, so that this hierarchy is in fact 
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integrable. This follows from taking time derivatives of (|5.10|) and using (|5.21|) to get 



d 



n-l 



dtr 



£=m—k 

which leads to the Zakharov-Shabat equations 

dQ m dQ 



n+k 
l=m+l 



0t r 



0t r 



+ 



Q m , Q n 







(5.26) 



(5.27) 



The matrix model actually defines a reduced, generalized KP hierarchy, because the flow 
equations are to be supplemented with the constraints imposed by the string equations. 

The knowledge of the solutions Q to ( |5.22| ) completely determines the partition func- 
tion of the fermionic matrix model. To see this, we consider the (n, n — 1) matrix element 
of the flow equation ( |5.24|) , which using (|5.13|) and (|5.15|) leads to the differential equation 

d\ogR n 



dt k 



Q 



Q 



n— l,n— 1 



(5.2* 



We can use ( \5.7\ ) and ( p.28| ) to obtain the flow equation for the partition function 

d log Z N 



dU 



Tr (iV ) (Q k ) 



(5.29) 



where we have defined the iV dimensional trace 

N-l 



TW°) = £[o] 



(5.30) 



n=0 



Using the KP equation ( p.22j ) and the Jacobi property ( |5.11| ) of the Q operator, we have 



g[Q]n 

dt k 



Q 



n+l,n 



Q 



n,n—l 



which, on using ( |5.28[ ) with k = 1, leads to 

d 2 hgZ N 



NN-l 



(5.31) 



(5.32) 



dtidt k 

This procedure can be iterated to determine any number of time derivatives of log in 
terms of the Q operators. Knowing Q therefore determines the free energy of the matrix 
model up to an overall integration constant. 



5.3 Fermionic r-Function 

It is clear from (|5.32|) that the partition function Z^ is a r-function for the discrete KP 
hierarchy (|5.22|) [j34|j . Let us introduce the Baker- Akhiezer functions 

m n [t-z\ = ^ ir ^y- (5.33) 
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and their conjugates 

e f v(z) A ( z ) 

y n [t ]Z ] = e -^ J ±P (5.34) 
which are bi-orthonormal in the standard line measure of the complex plane, 

£ dzV n [ t; z] ^ m [ t; z] = 5 nm (5.35) 

2 = 

The operator Q is a discrete Lax type operator which has eigenfunctions ^ n [t ; z] with 
eigenvalue z, 

Q* n [?;z] = z* n [?;z) (5.36) 

The function Z^[t ] is then the generating function for the Lax operator Q, in the sense 
that Q can be reconstructed from a knowledge of Z^. In fact, if we introduce the discrete 
r-function T n [t] for this hierarchy by the formula (|4.1|) ||15|| , then from fl3.7|) we see that 
the partition function of the fermionic matrix model is given by 

Z N [t}=r N [t] (5.37) 

Formally, the r-function is a section of the determinant line bundle over a Sato Grass- 
mannian associated with Riemann surfaces of the spectral parameters z ||34|| . It may be 



characterized as the unique function depending on the times tk and satisfying the given 
hierarchy of constrained differential equations. It is the generating function for solutions 
of the generalized KdV equations. Thus once the r-function is known, everything about 
the matrix model is likewise known. The usefulness of the relationship ( |5.37D is that there 
are a large number of identities satisfied by rjv[£ ] that can be used to characterize the 
partition function of the fermionic matrix model. 

The flow equations ( |5.22| ) or (|5.26| ) constitute the Lax representation of a discrete 



integrable hierarchy of differential equations of Toda type ||35|| . It is tempting therefore 
to characterize the fermionic partition function as the r-function of a certain reduction 
of the integrable Toda chain. However, as we will now demonstrate, this is not the case. 
By equating the constant terms on both sides of the flow equation ( |5.21| ) for k = 1 and 
using ( |3.26| ) and ( |3.27| ), we have 



^1 ^n+1,0 m= o k=m+l 

By iterating the relation ( |3.30|) using (|3.27| ) it follows that 



J] Rk (5.38) 



Pn, /n,0= (-l) n+1 (^n-l) (5.39) 

where we have defined R = 0. Substituting ( |5.39| ) into ( [5.38D we obtain 

^ = Pn+lj0 R n = Pn+h0 (l + (-1)" +1 Pnfl ln,o) (5.40) 
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Similarly, one can derive the evolution equation 

dl n ,0 



k-1,0 (l + (-l) n+1 Pn, /n,o) 



(5.41) 



The flow equations ( |5.40| ) and ( |5.41| ) are the first members of an integrable Hamiltonian 
system known as the Toeplitz chain hierarchy [p5fl . The general evolution equations may 
be similarly derived using the flow equations of the previous subsection and are given by 
the Hamilton equations of motion 



for the system of Hamiltonians 



01 



n,0 



dt k 



-l) n+1 Pnfllnfl) ^ (5-42) 



<9p 



'/ i.O 



#fc - — T LL (AT) 

and the symplectic structure 



A' 



1 <91ogZ 



fe > 1 



(5.43) 



(5.44) 



,o 



The Toeplitz chain is a particular reduction of the two-dimensional Toda lattice hierarchy, 



the latter of which characterizes generic matrix integrals |30| , |15|| . Thus the integrable 
hierarchy to which the adjoint fermion one-matrix model belongs is not the reduction 
of a one-dimensional (Toda) hierarchy, but rather of a two-dimensional one. This is of 
course anticipated from the lattice interpretation of the fermionic matrix model and its 
associated doubling of degrees of freedom due to the Penner matter states. 

However, the Toeplitz chain hierarchy can be viewed as a deformation of the standard 
Toda chain hierarchy. To see this, we use ( |5.28| ) for k = 1 and (|3.26|) for m = n to obtain 
the flow equation 

d log h n 



[Q]nn Qni^R 



n+1 



1 



On the other hand, by using ( |5.31|) for k = 1 and (|3.26|) for m 

d[Q]nn 



Qn Qn+lRn+l (R 



n+2 



n — 1 we have 

Qn-lRn(Rn+l ~ 1, 



(5.45) 



(5.46) 



The two flow equations ( |5.45| ) and ( 5.46|) may be combined together to give 
d 2 \ogh n dlogh n ( R n+ i d\ogh n+ i R n 91og/i n _i\ 



dt\ ' dti \R n+ i - 1 dtx R n -l 
By introducing time-dependent functions q n [t } through 

h n [t] = (-l) n e 2n e tl/£ e qn[t] 



(5.47) 



(5.48) 
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where e is a time-independent parameter, we can write ( |5.47| ) as 



d 2 q n 




-l-Qn 



dt\ \ dti J \ dt x J 1 + e 2 ei^+i-in 

pQn—qn-l 

—o (5-49) 

The differential equation ([5 .49 ) is the first member of an integrable hierarchy known as 
the relativistic Toda chain J36j. In the "non-relativistic" limit e — > 0, it reduces to the 
Toda chain hierarchy which characterizes generic Hermitian one-matrix models. Thus, 
the partition function of the adjoint fermion one-matrix model is a r-function of not the 
Toda chain hierarchy, but rather of its deformation to the relativistic Toda chain. The 
latter chain is itself a reduction of the two-dimensional Toda lattice. From this point of 
view, the fermionic matrix model may be thought of as a "deformation" of the Hermitian 
one-matrix model, the role of the deformation being played by the Penner interaction. 

The remaining differential equations of the relativistic Toda chain hierarchy are given 
by the Lax equations ( |5.26| ), by ( |5.28| ), and by exploiting the Jacobi properties ( p. lip 
and (|5.12|) of the operator Q. It is important to realize that, because the present model 
contains only a single set t of times, the ordinary Toda lattice does not itself appear in the 
hierarchy satisfied by the partition function. Thus the standard Toda lattice structure 
which underlies all integrable models (and in particular matrix integrals) only appears 
in a very subtle way through the reductions obtained above by eliminating one set of its 
time variables. This reduced structure leads to a certain degeneracy in the r-function as 
compared to the usual two-dimensional Toda lattice structure. Of course, the relations to 
integrable hierarchies described in this subsection are only 'kinematical', as they merely 
rely on the very basic recursion properties satisfied by the orthogonal polynomials. To 
incorporate the 'dynamical' aspects, and in particular the effects of the Penner potential, 
we need to impose the constraints implied by the string equations. Indeed, there are 
many r-function solutions of the above hierarchies, and the string equations select the 
one appropriate to describe the nonperturbative dynamics of the adjoint fermion one- 
matrix model. This is the topic of the next subsection. 



5.4 Virasoro Constraints 



We will now begin examining the constraints imposed on the fermionic r-function Zjq\t ] 
as dictated by the string equations. Proceeding as in ( |5.8|) for the operators PQ n+1 , 
n > — 1, we can express the string equation (|5.9|) as an infinite system of equations 



Tr 



(AO 



P + P f Q 2 - (N + 1) CT 1 + ]T kt k Q k 1 | Q 



ra+l 



fe>l 







n > -1 (5.50) 



The strategy now is to rewrite the equations of motion (|5.50| ) using the flow equations 
( |5.29| ) and ( p.32| ) . In this way, the string equations will be represented as the annihilation 
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of the partition function by a system of differential operators in the coupling constants of 
the potential V. In this subsection we will deal with the cases n > 0. 

Let us start with the n = string equation of (|5.50 ). Using the canonical commutation 
relations (15. 61) and the matrix elements 



h n [QP}nn = (z K( Z ) \ M*)) = ™K 



(5.51) 



we may compute the first trace in ( 5.50Q for n = as 



n=0 Z 



The second trace may be similarly computed by using 

h 



P f Q 1 



to get 



Tr {N) (p f Q" 1 



N(N - 1) 



-nh r 



(5.52) 



(5.53) 



(5.54) 



By substituting (p.52|) and (|5.54|) , and using ( |5.29|) , we can represent the n = constraint 
of the system ( |5.50| ) as the flow equation 



J2kt k ^-N 2 \Z N [t]=0 



(5.55) 



Next, let us consider the n = 1 constraint of (|5.50| ). For the first trace, we use ( j3~5l 
and ( |3.6| ) to compute the matrix elements 



h. 



Q 2 P 



z 2 &Jz) 



An(z) 

(nz n+1 - (n - l)p n ,n-iz n + O I A n (z)) 

np n+1 , n - (n - l)p„, n -i h n 



(5.56) 



Using ( |3.16| ) and the canonical commutator ( |5.6| ) we thereby arrive after a little algebra 
at 



N-l 



Tr (7V) P Q 2 = 2Tr (JV) (Q) + £ (n[Q] nn + Tr (n) (Q) = (N + 1) Tr (iV) (Q) (5.57) 



n=l 



Since 



= 0, upon substituting ( |5.57|) into (|5.50|) we find that the terms linear in Q 

n, 

cancel out and applying the flow equation ( |5.29|) we see that the n = 1 constraint can be 
written as 



k>l 







Of 



k+1 



z N [t] = o 



{5.5i 
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Now we move on to the n = 2 string equation of ( |5.50| ). Again using (|3.5|) and ( ^6| 
we may compute 



h, 



Q 3 P 



z 3 &Jz) 



(n Z n+2 _ {n _ l)pnn _ x z n+i _ {n _ 2) Pntn _ 2Z " + O (z^ 1 ) | A n {z)) 



n(p n+ 2,n + Pn+2,n+lPn+l,n) ~ (u - 2)p n ^ 2 ~ (n - l)p n+ i jn p n>n _i 



h, 



(5.59) 



n-l 

Pn,n-2 = X]([Q]/ 
k=l 



(5.60) 



Using ( |3.16| ) we can iterate the relation ( |3.14| ) for k = n — 2 to get 

U-./,;-l + [Q]kk Tr(fe)(Q) ; 

Substituting ( |5.60| ) and ( |3.16j ) into ( |5.59| ) and using (|5.6|), we arrive after some algebra at 
Tr (7V) (PQ 3 ) = 3Tr w (Q 2 ) + (iV-l)[Q]^_ 1 + (iV+l)(Tr (i v ) (Q)) 2 

N-l 

+ {2N-3) J2{[Q]n,n-l-[Q}nnTr {n) {Q)] 



(5.61) 



n=l 



The sums in (|5.61|) can be simplified by using the Jacobi property (|5.11|) of the operator 
Q to get 

V ' A 2 



Tr w (Q 2 ) = [Q] 



N-l N-l 
N,N—1 2 ^ f [Q]n,n— 1 

+ E ([Q] 

n=l n=0 



(5.62) 



and in this way we arrive at an expression for the first trace in the n = 2 equation of 
Tr (7V) (P Q 3 ) = (# + §) TV) (Q 2 ) + \ [Q]n,n-i + \ (Tr (iV) (Q)) 2 (5.63) 



For the second trace in ( |5.50| ), we use ( |3~3| ) to compute the matrix elements 



P Q 



z 3 K(z: 



($ n (z) I -n z~^ + O (^("- 3 )) + l nA z) 



L i (*$ n (z) 



Using ( |3.25|) we obtain 
where we have defined 



P T Q =S n Q n {l-R n+l ) 



C _ ^,1 



(5.64) 

(5.65) 
(5.66) 



n,0 



By equating the coefficients of the z~ x terms on both sides of the three-term recursion 
relation ([3.32 ) for the A polynomials, we can write down an iterative equation for the 
coefficients ( |5.66|) . 

S n = RnQn-l + S n+ \ — Q n (5.67) 
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which has solution 



ra-1 



&n Qn—l 



(5.68) 



fc=0 



Substituting ( p.68j ) into fl5.65| ) and using ( |3.26| ) for m = n — 1 and m = n, we therefore 
find 

"P^Ql = [Q]nnTr (n) (Q) - [Q]„,„_i (5.69) 

L J nn 

Using the identity ( [5.62 ) we then arrive at an expression for the second trace in the n = 2 
equation of ( |5.50| ), 



Tr 



(AT) 



P f Q) = \ [Q]jv,iv-i + \ (IWQ))' - 1 Tt (n) (Q 



(5.70) 



Upon substitution of (|5.63| ) and (|5.70| ) into ( |5.50| ) for n = 2, we see that the terms 
quadratic in Q cancel out. The remaining terms can be simplified using the flow equations 
( |5.29|) and ( |5.32| ), which leads to the n = 2 constraint equation 



Of 



fc+2 



dt\ 



z N [t 







(5.71) 



This procedure can be easily generalized to higher orders n > 3. However, it is possible 
to determine the general relation from the first three equations (|5.55|) , ( ^.58| ) and ( ^.71 ) 
by using the fact that the commutator bracket of the pertinent differential operators 
must also annihilate the partition function and thereby using commutators of the above 
operators to generate the higher order ones. The final result can be expressed in a more 
familiar form by introducing the zeroth time to through the flow equation 

dZ N 



dt. 



NZ 



N 



(5.72) 



for the partition function. 



In this way we find that the string equations ( |5.50| ) for n > can be written as the 
set of discrete Virasoro constraints 







n > 



where the second order linear differential operators 

d 



k>l 



dt 



k+n k=0 



" 



dtu dt 



2N 



n—k 



d_ 

dt„ 



satisfy the Witt algebra 



(n — m) L r 



(5.73) 



(5.74) 



(5.75) 



for n,m G Z + . These Virasoro constraints are the usual ones for the fermionic matrix 
model ||, Q and they are identical to those of the Hermitian Penner matrix model [37| 
as defined in (Of) with a = —2. They represent the full system of equations of motion 
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of the matrix model and can be equivalently derived by demanding the invariance of 
the matrix integral under arbitrary changes of the matrix variables || ||. The first set 
of operators in ( |5.74 ) comes from the variation of the potential V in the action, while 



the second one comes from the Jacobian of the change of matrix integration measure. 
The first two terms in ( |5.74|) therefore coincide with the standard Virasoro generators of 



generic Hermitian one- matrix models [|30|. The last operator in ( |5.74| ) comes from the 
variation of the logarithmic Penner interaction. 

It is a remarkable fact that the loop equations of the adjoint fermion one-matrix model 
( |1 . 1| ) , the Hermitian Penner one- matrix model (( |1.3| ) with a = —2), and the unitary 
Penner matrix model ( |1.4|) are all equivalent. However, the equivalence between the 
Hermitian model and the other two matrix models only holds order by order in the 
expansions of the matrix integrals. Although the models are equivalent at N = oo, beyond 
leading order in the large N expansion the loop equations should be solved with different 
boundary conditions and the solutions are different. The fermionic and unitary equations 
should be solved with boundary conditions appropriate to a perturbative Gaussian limit, 
in order that the models admit interpretations in terms of random surfaces, while the 
Hermitian ones should be solved with boundary conditions appropriate to recover the 



Penner model of the discretized moduli space of Riemann surfaces |10|, jlgl- The distinction 
of the fermionic and unitary matrix integrals from the Hermitian one is evident from the 
structure of the finite N correlators (see appendix A). In fact, it is only in the large N 
limit that a full Virasoro symmetry is realized in these models, as well as an infinite 
hierarchy of generalized KP differential equations. This makes the large N limit a very 
natural ingredient of the fermionic matrix model, as is also apparent from its fat-graph 
interpretation of section 2.2. 



5.5 Origin of the Painleve Expansion 



In the previous subsection we have shown that the Schwinger- Dyson equations of the 
fermionic matrix model are generated by the positive Borel subalgebra of a Virasoro 
algebra of vanishing central charge. There is a nice algebraic way to characterize the 
effect of the Penner interaction potential in the Virasoro generators ( |5.74j ). Let us write 
them as L n = + T n , where are the standard Virasoro generators of a Hermitian 
one-matrix model with potential V, and 





T„ = -2N 



dt r 



(5.76) 



are the generators of translations t n *—>■ t n — 2Ne in coupling constant space. Together 
these operators generate the symmetry algebra 



T T 

± n i ± rr 



Tti Tt 







In 



m)Ll 



+m 
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T T u ~ 

± n -i lj m. 



-2NnT n , 



(5.77) 



The commutation relations ( 5.77 ) characterize a master symmetry of the generalized KdV 
hierarchy ]33J. We see then that the Virasoro constraints of the fermionic matrix model 
constitute a deformation of those for Hermitian one-matrix models by the master sym- 
metry algebra generators of the KdV flow equations. Note that the particular symmetry 
generated by the operators ( |5.76|) is restricted to translations in units of a discrete "lat- 
tice spacing" —2N, indicating what sort of reduction of this master symmetry is taken. 
These Virasoro constraints describe the desired reduction of the fermionic r-function, and 
they lead immediately to the desired geometrical interpretation of the Painleve expansion 
of the adjoint fermion one-matrix model. The Painleve equations may be characterized 
as following from a reduction of the r-function of the integrable Toeplitz (or relativistic 
Toda) chain hierarchy, subject to the Virasoro symmetry obtained in the previous subsec- 
tion. They yield a stability condition L n = on the points of the model Grassmannian, 
which are associated with the Baker- Akhiezer functions ( |5.33| ), that select a particular 
class of transcendental solutions to the KP equations. 

It is a standard fact that a set of operators of the form ( 5.74j ) can be embedded into 
a Virasoro algebra of non-vanishing central charge. However, an important aspect of the 
present approach to the fermionic matrix model is that there are in fact extra dynamical 
constraints imposed on the partition function which can be attributed to the invertibility 
of the Lax operator Q. For instance, the constraint ( |5. 50 ) makes perfect sense for n = —1, 
and the calculation of the trace of the operator P Q~ 2 proceeds in an analogous manner 
to that of (|5.56|) with the result 



TV) (P f Q- 2 ) = ~(N - 1) Tr (A0 (Q- 1 ) 



(5.78) 



Since [P] nn = 0, we see that on substitution of ( |5.78| ) into (|5.50|) for n = — 1 we produce 
an extra non-vanishing term — 2N Tr^^Q -1 ). This trace can be written in an operator 
form as follows. From ( |3.31| ) with m = n and ( |5.45| ) we have 



Q 



Rn-l) R 



1-71+1 



1 



'dlog h n 



(5.79) 



This quantity is independent of the first time t\, since on using the flow equation fl5.24|) for 
k = 1 and the Jacobi properties of the operators Q and Q we find d[Q _1 ] nn /d£i = 0. This 
suggests defining a negative time t_ x such that Tr^^Q -1 ) is represented as an operator 
acting on the partition function as 



DZ 



N 



d DZ N 
dh Dt-i 



N-l 
n=l 'Wn+l 



1 (t„+iT„_i - 7-2) (r n+2 T n - T, 



2 

n+1 



dr n+ i 



T , din 

ln + l at! 



(5.80) 
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where we have used ( |4.1|) . Then the constraint (|5.50|) for 



n 



T, kt kT^— + Nt 1 -2N-^- ) Z N [t 

Otfc-l JJt-i 



-1 can be written as 
= (5.81) 



> k>2 



In fact, there are infinitely many negatively moded constraints such as this, associated 
with higher powers Q~ m , m > 1. However, the operators such as that in ( 5.81j) which are 
thereby generated do not form a closed algebra among themselves or with the operators 
(|5.74j) . Indeed, this set of constraints arises from the previously mentioned reduction from 
a two-dimensional lattice (containing both negative and positive time parameters) to a 
(deformed) chain by eliminating all negative couplings. As usual, such reductions break 
the full Virasoro symmetry of the original integrable model and only the positive Borel 
subalgebra survives here as a symmetry of the matrix model. The negatively moded 
constraints are therefore merely an artifact of the lattice reduction or equivalently the 
chain deformation. In practical terms, all observables in the fermionic model are described 
by observables conjugate to positive times in the unitary matrix model. Observables 
conjugate to negative times exist but do not appear to play a role in the equivalence. 

We can understand all of these matters more clearly by rewriting the Virasoro con- 
straints Q5.73 ) in terms of the Baker- Akhiezer functions (|5.33|) . For this, we note first of 
all that ( |5.21 ) and ( 5.28|) imply that they obey the flow equations 



d* n [t; z] 
dtk 



Q*- Q 



^n[t; z] 



(5.82) 



Now let us compute the action of the scale transformation generator z in eigenvalue 
space on the Baker- Akhiezer functions. We have 

99 n [t;z] i .V 3 



N \ ef y W 

zV'(z))^ n [t;z] + ^ ir PQ^ r 

A / 7. 2 



(5.83) 



dz V 2 & / z - 

To evaluate the last term in ( 5.83|) , we have to take into account the triangularity of the 
operator P Q. From P = P_ and the Jacobi property (|5.11|) , we have 



PQ = (PQ)_ + £[P Q] nn E n>n 



(5.84) 



n>0 



We can therefore evaluate the action of P Q on & n (z) by using the string equation ( |5.9|) 
multiplied on the right by Q, but eliminating the pure upper triangular part of both sides 
of the equation. Using ( |5.53| ) and the fact that 



t ^ i 



o 



(5.85) 



which can be derived similarly to the other matrix elements of this section, we can write 



(I5783D as 



d^ n [t; z] 

dz 



(2n + N-i)V n [t; z] 



N 



(q v'(Q)) - (q V(Q))_ - [Q v^'(q) 



tf n [f;z] (5.86) 
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The last term in (|5.86|) can be rewritten using the flow equation ( |5.82|) , leading to the 
eigenvalue equation 

L [t;z]^ n [t;z] = -n^ n [t;z] (5.87) 
where we have defined the Virasoro operator 



d d N — 1 



(5i 



k>l 



The eigenvalue equation (|5.87|) shows that the lowest Virasoro constraint is indeed a 
true symmetry of the system, in that the operators Q and Lo[£;z] are simultaneously 
diagonalizable in the basis of bi-orthonormal Baker-Akhiezer functions. The same is 
true of all higher Virasoro constraints. This shows precisely what sort of constraints 



the auxiliary eigenvalue problem (|5.36|) must possess in order to correctly reconstruct the 
fermionic r-function. However, an analogous computation using the translation generator 



on eigenvalue space leads to the result 



v fc>i 



d 



dt 



k-1 



d_ 

dz 



N + l 

2z 



V n [t; z] = -{n + N)R n #„_i[f; z\ 



(5.89) 



which shows that there is no corresponding L_i operator which commutes with the Lax 
operator Q. This means that there is no translational symmetry in the system, and 
the leading equations of motion at large N are determined instead by the generator Lq of 
scale transformations of the system. These equations can be written in a form comparable 
to those of sections 3 and 4 by writing a consistency condition between the eigenvalue 
equations ( |5.36| ) and (|5.87|) in the form 



(Lo[t; z]+n)(Q-z)V n [t; z] = 



(5.90) 



Expanding ( 5.90|) using ( |5.87| ), (|5.88|) and completeness of the Baker-Akhiezer functions 
then leads to the system of equations 



(1 + n-m) [Q) nm + E kt > 



d[Q], 



k>l 







(5.91) 



In particular, setting n = m in ( |5.91[ ) and using the flow equation ( |5.31| ), we have 



[Q]t 



E kt k ( [Q 



k>l 



n,n—l 



n+l,n 



(5.92) 



What this means is that it is the operator P Q (along with its conjugate Q -1 ) 
which plays the fundamental role in the dynamics of the adjoint fermion one-matrix 
model. Thus it is the scale invariance of the system which leads to the fundamental 
double scaling equations in the large iV limit. The translational symmetry generated 
by the operator P itself (and its conjugate P Q™ 2 ) is an extra constraint on the theory 
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which plays no immediate role in the continuum limit of the matrix model. By considering 
the matrix elements {& n (z)\z [zA m (z)]), it is possible to infer the commutation relation 
Q , P Q" 1 = Q. Using the canonical commutator (|5.6|) , it follows that the fundamental 
symmetry operator of the system is given by 

n=i(PQ-P t Q- 1 ) (5.93) 



which together with the Lax operator Q obeys the non-canonical commutation relation 

n,Ql=Q (5.94) 



The operator expression ( |5.94p defines the appropriate string equations of the fermionic 



matrix model. The upper and lower triangular parts of the scaling operator (|5.93|) may 
be computed by using ( |5~9D , (|5.51|) , (|5.84|) and (|5.85|) to get 



n + = \ J2( n - N ) E n, n + l E fc ^Q+ 

A n>0 A k>l 

U = - l -Y,kt k Qt (5.95) 

1 k>\ 



Substituting ( |5.95| ) into the string equation ( |5.94j ) and using the KP equation ( J5.22Q , we 



arrive at (|5.91| ) 



It is the equation (|5.92| ) which produces the novel Painleve expansions of the fermionic 
matrix model. For instance, in the case of the quadratic potential studied in section 4, 
it is implicit in the manipulations carried out in section 4.2. In fact, the structural form 
of ( |5.91| ) is identical to the so-called "automodel" constraints which arise in the usual 
Hermitian one-matrix models [^] . On the other hand, in the continuum limit Q becomes 
a differential operator of a certain finite degree and the string equations (|5.94j) can be 
translated into an equation for pseudo-differential operators which satisfy a generalized 
KdV hierarchy of flow equations [|]. General solutions of string equations such as ( p.94| ) 



have been studied before, in the case that the double scaling limit of Q is described 
by a Schrodinger operator, with stable, pole-free solutions for the corresponding string 
susceptibility |38| . We see therefore that the string equations of the fermionic one- matrix 
model lead to a Painleve differential equation which has the same structure as that in the 
usual Hermitian one-matrix models, except that now the automodel form of the equations 
produce an alternating genus expansion. Moreover, the structural form (|5.94|) of the string 
equations allow us to restrict to positive values of the double-scaling variable £ and thereby 
obtain pole-free solutions for the free energy of the fermionic matrix model. 

In the present case, the details of the continuum limit of the operator Q are somewhat 
involved. However, the discrete equations Q5.91D , Q5.92| ) and ( |5.94j ) completely characterize 



the double scaled equations for the partition function in the large N limit and the ensuing 
topological expansion. They serve as the starting point for a characterization of the 
differential hierarchies of the fermionic one-matrix model. From the generic structure of 
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the flow equations obtained in this section, we can expect that the continuum equations 
satisfied by the partition function of "fermionic quantum gravity" are closely related to the 
usual KdV flow structure of two-dimensional quantum gravity in terms of Gelfand-Dikii 
differential polynomials 0. |j . 
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Appendix A Spectral Correlation Functions 

In this appendix we will briefly describe some properties of correlators of the fermionic 
one-matrix model within the formalism of this paper. The symmetries of the matrix 
integral ( |1.1|) restrict its observables to those which are invariant functions of ipip 0. 
Connected correlation functions may be generated by taking derivatives of the free energy 
with respect to the coupling constants of the potential V, 

F , conn J J 

where the normalized fermionic correlation functions are defined by 

(/ (^)) F = / # d$ f eF*v&» (A.2) 

Gr(7V) c 

and in ( |A.1[ ) it is understood that, if necessary, a set of auxiliary coupling constants 
Qk are introduced into the potential V and then set to zero after differentiation. Given 
the equivalence of the generating function logZ/v for the connected correlators of the 
fermionic and unitary matrix models, we therefore also have complete equivalence of their 
observables, 

/ l \ 

pj \ 






(A.3) 

F,conn 1 ' U,conn 

where the normalized unitary correlation functions are defined by 

(f(U)) v = ^L J [dU] /(CO e^(^)- 1 ^) (A.4) 

N U(N) 
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The correlation functions in (|A.3|) may be evaluated as in ( p.29|) -( p.33|) . For example, 
one may readily compute 



tr 



e ip9 trUg 



NZ 



N 



[dU ] trU p e NtlV ^ 



SU(N) 



(A.5) 



d=N 2 -p 



The restriction in (|A.5|) of the character expansion of e NtTV i u o) to Young tableaux with 
A^ 2 — p filled boxes implies that single trace correlators are non- vanishing only for p < N 2 . 
The finiteness of the set of non-vanishing correlators is of course natural as a consequence 
of the anticommuting property of the matrices ip and ip, but it is a dramatic result in 
the unitary matrix model. Thus not only does the Penner interaction term in the unitary 
matrix model (|2.28|) capture the finite nature of the perturbation expansion of Z N , but 



it also reduces the correlation functions appropriately to reproduce the correct properties 
of the fermionic correlators at finite N. Note that, geometrically, the operator tr(if)ip) p 
inserts, on the dual triangulated surface, a hole with p boundary lengths. The lattice 
expansion of such operators thereby generates fermionic ribbon graphs which are dual to 
tesselations of Riemann surfaces with a given number of boundaries. Again there are only 
finitely many such fat-graphs, because the maximum number of boundaries that a given 
discretization can have is A^ 2 . 

The complete set of observables of the fermionic matrix model can be generated by 
the joint probability distributions 

(N-n)\ 



Z n ) 



N\ 



11 tr 5 (zk-Hf) 



(A.6) 



\fc=l / f 

where Zk are points in the complex plane and 1 < n < N. When the correlation function 
in ( |A.6| ) is mapped to the unitary matrix model as described above, the points Zk can be 
interpreted as eigenvalues of unitary matrices. By diagonalizing the unitary matrices in 
(|A.4|) , using the identities 



A(zi, ...,z N ) 



det 

i,3 



A(zi, ...,z N ) 



det 



Aj-i(zi) 



(A.7) 



which follow from (|2.6|) and ( |3.5| ), and by using Q3.7| ), it is straightforward to derive the 
determinant representation 

(N-n)\ /" \ (N — n)\ 



Pn\Z\i ■ ■ ■ i Z n ) 



Y[tr5(z k -U) 



\k=l 



det 



where )C(z, z') is the spectral kernel which is defined in terms of the orthogonal polynomials 



as 



1C(z, z ) 



\zz' 



JV+1 



n=0 



K. 



(A.9) 



We see therefore that the problem of evaluating correlation functions of the fermionic 
matrix model reduces to that of determining the spectral kernel ( |A.9| ). It is possible 
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to express it in a much simpler form by deriving the appropriate generalization of the 
Christoffel-Darboux formula [|l|. For this, we first need to derive a "mixed" recursion 
relation between the $ and A polynomials. Consider the bi-orthogonality relations 



$ n+ i(Y) -z$ n {z) 
z n A n (z) 



$n+l(*) 



$n0) 



-(m-1) 







A n (z) 







(A.10) 



which are valid for 1 < m < n. From ( |A.10| ) it follows that the two polynomials § n+ \(z) — 
z Q n (z) and z n A n (z) of degree n are equal up to some constant c n . Equating the constant 
terms of these polynomials gives c n = p n +i,o, and we arrive at the mixed three-term 
recursion relation 

$ n+ l0) = Z $ n (z) + Pn+1,0 Z n A n (z) (A.ll) 

In an analogous way, it is possible to derive the recursion relation 



A„+iO) = - A n (z) + l n+1 o \ § n (z) 



(A.12) 



We now multiply ( |A.11| ) through by A n (z')/h n and ( |A.12| ) with the index shift n — » 
n — 1 and z — > z' through by z' $> n (z)/h n , subtract the resulting two equations, and then 
sum over n = 1, . . . , N — 1. This yields 



N ~ x $ n (z) A n (z') 



$ N (z) An-Az') $i(z) 



n=l 



h. 



N-l 1 
n=l n k 



hi 

Pn+lfi Z n A n ( Z ) An(z') - -^T * n {z) ^n-l(z') 



(A.13) 



We can iterate the recursion equations ( |A.11|) and ( |A.12| ) to get 



A fc (z) = l + *£ij,o$j-i(z) 



fc-i 



y/fc 



Substituting ( |A.14 ) into ( A.13 ) and comparing the resulting equation with itself under 
the interchange of arguments z <-> z', we arrive after some algebra at 



^ _ z >) J] ?»(*) A ^ ') - ^A^^^W-^A^!^)^^') 



n=0 



/ijV-i z 



/N-l 



+ — (2 - 2 ) 

h 



(A.15) 



The expression ( A.15 ) is valid for any pair of complex numbers z ^ z' . We can take the 
z — > z' limit of (|A.15) using l'Hospital's rule to get 



N ^ $ n (z)A n (z) _ ^ JV ~ 1 A N ^i(z) <&' N (z) - &n(z) (z^ 1 A n ^(z))' i- ho 



E 

n=0 



(A.16) 



53 



The identities QA.15 ) and (|A.16 ) are the fermionic analogs of the Christoffel-Darboux 
formula for the usual Hermitian one-matrix model orthogonal polynomials. 

The Christoffel-Darboux formula allows us to express the spectral kernel (|A.9|) in 
terms of orthogonal polynomials with indices n close to the matrix dimension N. Using 
the mixed recursion relation (|A.ll ), we obtain 

e f (v(z)+V(z>)) 



fC(z, z) 



ZZ' 



N+l 



Z 



l- N+l 



h N -iPN,o {z - z') 



Z <$>N-l{z) ^n(z') - Z ^N-l(z') ® N (z) 



1-hr 



(A.17) 



It follows that the fermionic matrix model is completely determined by the single set of $ 
polynomials, as expected because of the one-matrix nature of the model and the heuristic 
identification U ~ ipip. In particular, for the spectral density we find 



p(z) = -jf]C{z,z 

e NV(z) ( 



Nz N {h N _ lPNfi z N 
I -h ] 



z (<$> N (z) $' 



JV-lv 



+ 



h 



- & N (z) $N-l(z)) + $ N (z) ^N-i(z) 



(A.18) 



The formula ( |A. 18| ) is particularly useful for determining the spectral density, and 
hence all correlators, of the fermionic matrix model in the large N limit. In that case, 
factorization and symmetry imply that all connected correlation functions vanish and the 
large N limit of the model is completely characterized by the set of correlators 



-tr 



dz p(z) z p 



(A.19) 



where the integral goes over the support of the function p(z) in the complex plane. Note 
that the spectral density is defined a priori in the adjoint fermion matrix model by the 
formula ( |A.(j| ) [TT| , |J. In the present formalism it has a natural interpretation as the 
probability density of eigenvalues in the corresponding unitary one-matrix model. As 
such, it is generally supported on the unit circle. However, as discussed in the previous 
subsection, the restriction of the integration contour to the unit circle in the large N limit 
is not necessary, and the spectral density can be supported generically in the complex 
plane. These facts explain the general properties of the spectral densities of fermionic 
matrix models l 



11 



], for instance how the pole generated by the Penner interaction 
requires the support contour of p(z) to be adjusted so as to avoid the origin of the complex 
plane and why the support endpoints are complex-valued. 
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Appendix B Solution of the Gaussian Model 



In this appendix we will illustrate how the formalism of section 3 works in a specific 
example. We shall consider the example of the Gaussian potential 



V(z) = z 



(b.i; 



for which everything can be obtained explicitly. In this case the recursion relations of 
section 3.2 are easily solved to give the coefficients 



Rn 
Pn 
Qn 
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n + N 
n 



N 



1 + 



n+ 1 
N 



and the normalization constants are 

h„ = 2ni N N 



nl 



n + N)\ 

Substituting (|B.2|) into the three-term recursion relation (|3.29| ), we obtain 



:b.2) 



(B.3) 



n-l + N\ (n-l)z 



N 



N 



(B.4) 



To solve the recurrence relation (|B.4|) , we introduce the generating function which is 
defined as the formal power series 



(B.5) 



n=Q 



1 



in a variable s, with the boundary condition En{z; 0) = 1. The recursion relation ( p.4j ) is 
then equivalent to a first order inhomogeneous linear differential equation for the function 

" Ar ' r 2 

& N ^ ~ SZ ) ds 

Integrating ( |B.6|) we find that the generating function (p.5|) is given by 

fc+i 



(l-sz)- + s(z + l 



SZ 

~N 



„ i ~ (N + k-iy. ( i \ k 

^ N {z;s) = -}^ 



k=0 



(N-l)\ 



1 — SZ 



(B.6) 



(B.7) 



Expanding the function ( |B.7|) as a power series in s and equating its coefficient of s n with 
that of the definition ( |B.5| ), we arrive at an explicit form for the polynomials <& n (z), 



(B.8) 



« fn\ (N + k- 1)! / 1 \ k n „ k 



k=0 



55 



which can be expressed as 



$ n ( z ) = N N z N+n U[N,n + 1 + N;Nz] (B.9) 

where 

7r / iFi[a, c;x] x 1_c iFi[a — c + 1, 2 — c; x] 
[a ' c;;rJ = ^r(a-c+l)r(c) + r(a)r(2 - c) 

is a confluent hypergeometric function. The partition function and observables of the 
Gaussian fermionic one-matrix model are therefore completely determined by a system of 
confluent hypergeometric polynomials. 



(B.10) 
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